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FOREWORD 


Recent  progress  in  instrumentation  and  computing  methods  suggested  that 
satellite  geodesy  ought  to  account  for  even  such  minute  effects  as  the  surface 
displacement  caused  by  the  tides  of  the  solid  earth.  To  reflect  the  earth’s 
deformation  due  to  the  gravitational  action  of  moon  and  sun  in  our  satellite  work 
at  the  Naval  Surface  Weapons  Center  (NSWC),  a suitable  physical  model  was 
established  for  the  tidal  surface.  Based  on  this  model,  a computer  algorithm  was 
written  which  specifies  the  tidal  perturbation,  in  position,  for  any  geodetic  station 
on  the  earth’s  surface  and,  in  particular,  for  the  satellite  observing  stations.  The 
numerical  trial  runs  resulted  in  surface  displacements  which  were  indeed  significant, 
and  the  just  mentioned  algorithm  was  subsequently  included  into  the  CELEST 
system  of  computer  programs. 

The  work  documented  here  was  done  in  the  Astronautics  and  Geodesy  Division 
and  was  funded  as  part  of  the  CELEST  development  effort.  Mrs.  Louise  Gordon  of 
the  Computer  Program  Division  coded  the  algorithm  and  made  the  required 
computer  runs. 

TTiis  report  was  reviewed  by  Richard  J.  Anderle,  Head,  Astronautics  and 
Geodesy  Division. 
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INTRODUCTION 

CELEST*  is  a new  system  of  computer  programs  to  calculate  precise 
ephemerides  and  associated  diagnostic  data  for  artificial  earth  satellites.  Main  data 
source  are  the  observations  conducted  at  a worldwide  network  of  geodetic  stations 
of  the  Doppler  frequency  shift  on  radio  signals  transmitted  by  the  orbiting  satellites. 

CELEST  goes  to  a great  length  to  accomplish  its  accuracy  goals.  For  example, 
dynamical  accuracy  is  assured  by  a force  model  for  the  satellite  motion,  which 
accounts  not  only  for  the  customary  causes  of  orbit  perturbations  such  as  a very 
detailed  earth  gravitation,  air  drag,  and  light  pressure  but,  in  addition,  specifies 
several  rather  minute  perturbing  accelerations  like  those  caused  by  the  gravitational 
effects  of  the  tidal  deformations  experienced  by  the  ocean,^  the  atmosphere,'*’*  and 
the  solid  earth.®  For  the  sake  of  geometrical  accuracy,  there  is  a flexible  refraction 
subroutine  in  that  part  of  the  physical  model  which  formulates  the  radio-wave 
propagation  from  the  orbit  down  to  the  individual  observing  station.  State-of-the-art 
values  are  used  for  the  observing  station  coordinates,  and  just  like  the  dynamical 
part  of  CELEST  was  augmented  by  certain  tidal  terms,  there  was  recently  a tidal 
correction  attached  to  the  range  and  range-rate  data.  This  correction  introduces  the 
displacement  due  to  the  lunar  and  solar  tides  of  the  solid  earth  of  the  surface 
under  the  observing  stations.  The  nature  of  this  correction  shall  be  discussed. 


DETAILS  OF  THE  PHYSICAL  MODEL 


LIMITING  ASSUMPTIONS 

The  computer  algorithm  below  was  derived  under  the  following  assumptions. 

The  Love  numbers  are  constants.  Only  those  Love  numbers  that  regulate  the 
surface  deformation  appear  in  the  algorithm.  They  shall  be  named  “geometrical” 

Love  numbers  and  shall  be  symbolized  by  h^  to  distinguish  them  from  the  j 

“dynamical”  Love  numbers  k.  that  will  appear  closer  to  the  end  of  this  report.  1 

Our  algorithm  is  based  on  the  customary  linear  relationships  amongst  the 
harmonic  constituents  of  the  tide  potential  and  the  corresponding  constituents  of  the 
surface  deformation.  Each  of  the  constants  appearing  in  these  relationships  is  one  of 
the  geometrical  Love  numbers.  The  just  mentioned  relationships  themselves  are  being 
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accepted  without  any  detailed  scrutiny  of  the  relevant  geophysical  elasticity  theory 
(Reference  7,  Chapter  VII;  Reference  8,  Chapter  3,  Paragraph  III;  Reference  9). 

Tlie  earth  is  approximated  by  a sphere  of  radius  R.  R is  the  semiinajor  axis  of 
some  suitable  reference  ellipsoid,  the  particular  choice  of  which  is  uncritical.  All 
terms  are  neglected  which  are  of  the  order  of  the  square  of  the  reference  ellipsoid 
eccentricity  or  of  higher  order.  This  assumption  is  valid  because  the  surlace 
distortion  produced  by  our  theory  will  not  need  to  be  known  to  better  than  a 
centimeter,  which  implies  that  no  more  than  two  significant  digits  will  be  needed. 

Only  the  radial  component  of  the  tidal  deformation  of  the  land  surface  is 
significant.  This  can  be  rationalized  from  the  fact  that  only  a few  terms  will  be 
considered  in  the  surface  harmonic  expansion  of  the  surface  displacement.  Of  these, 
only  one  term  will  actually  appear  in  the  algorithm  because  it  alone  is  sufficient  to 
produce  centimeter  accuracy.  These  terms  (and  especially  the  one  actually  used)  have 
a very  long  wavelength  (thousands  of  kilometers)  and  a rather  small  amplitude  (a 
few  centimeters  at  most),  causing  the  surface  tilt  to  be  negligible. 

The  study,  of  which  the  present  report  is  a summary,  was  strictly  limited  by 
the  requirement  to  produce  a computer  algorithm  for  the  displacement  of  stations 
on  the  surface  by  the  tidal  mechanism.  This  objective  was  rather  single-mindedly 
pursued.  In  particular,  no  attempt  was  made  to  formulate  the  tidal  surface  tilt  for 
possible  use  with  certain  classes  of  observational  data  sensitive  to  the  tilt  angle  (for 
example,  astronom  eal  latitude  or  zenith  angle).  Neither  was  it  felt  advisable  to 
introduce  ocean  tide  loading  at  the  present  time. 

The  sun  and  moon  positions  are  specified  by  the  contents  of  certain  ephemeris 
computer  tapes  and  by  an  interpolation  routine  for  the  tapes.  These  ephemerides  are 
readily  available  because  they  are  frequently  used  elsewhere  in  CELEST. 


To  facilitate  those  derivations  which  concern  the  tide  potential  and  the  tidal 
surface,  it  was  at  first  assumed  that  the  tide  is  an  equilibrium  tide.*  As  explained  in 
detail  below,  the  missing  tide  lag  was  introduced  later  by  interpreting  the 

equilibrium  tide  potential  (and  also  the  resulting  surface  distortion)  as  being 
produced  by  a fictitious  moon  and  sun.  It  will  be  seen  how  a fairly  simple 

transformation  from  the  lunar  and  solar  ephemeris  coordinates  to  those  of  the 
fictitious  moon  and  sun  will  change  the  equilibrium  potential  field  and  the 
equilibrium  tidal  surface  into  a potential  and  a surface  which  contain  the  desired 
effects  of  tidal  lag. 

equilibrium  dde  is  a stationary  solution  of  the  mathematical  tide  problem.  It  is  different  from  the  real  tide  mainly 
because  it  wores  the  effects  of  the  inertia  and  viscosity  of  the  earth’s  masses  which  in  the  real  world  would  cause  a time 
W to  exist  betwee-  the  motion  of  sun  and  moon  in  the  sky  and  the  mass  redistribution  occurring  as  a consequence  within 
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EXPANSION  OF  THE  LUNAR  GRAVITY  POTENTIAL 


With  the  sole  exception  of  the  thermal  component*  of  the  air  tide,  ail  earth 
tides  are  in  essence  mass  relocations  within  the  solid  earth,  ocean,  and  atmosphere 
caused  by  lunar  and/or  solar  gravitation.  To  be  precise,  the  earth  as  a whole  rotates 
within  the  fields  of  lunar  and  solar  mass  attraction.  These  fields  impart  on  each 
terrestrial  particle  an  acceleration  which  is  a function  of  the  position  of  the  particle 
within  the  field.  Generally,  the  lunar  and  solar  force  fields  are  variable  across  the 
earth.  They  also  depend  on  time  in  accordance  with  tlie  motions  of  the  two 
celestial  bodies  relative  to  each  particle.  The  acceleration  of  the  particle  relative  to 
the  earth  as  a whole  in  turn  gives  rise  to  certain  geophysical  forces  which  balance 
it,  thus  enabling  the  particle  to  remain  attached  to  the  earth.  The  latter  forces 
constitute  the  tidal  stress  field  which  is  responsible  for  the  tidal  deformation  of  the 
earth  and,  in  particular,  for  the  distortion  of  the  earth’s  surface. 

First,  we  shall  discuss  the  lunar  field  in  the  earth’s  vicinity.  The  next  paragraph 
will  introduce  the  tide  potential  associated  with  the  lunar  gravitation.  The  TIDE 
POTENTIAL  AND  SURFACE  DEFORMATION  section  will  establish  the  surface 
deformation  (tidal  bulge).  Then,  the  tidal  lag  will  be  formulated  in  close  connection 
with  the  procedure  for  lunar  position  and,  finally,  the  solar  contribution  to  the  tide 
will  be  explained. 

Let  y^'\  i = 1,  2,  3,  be  an  earth-fixed  Cartesian  coordinate  frame,  i = 1 and 
i = 2 indicate  the  equator;  y^’^  and  y^^^  define  the  Greenwich  meridian.  COE  is 
the  center  of  the  earth  in  the  sense  of  CELEST  and  TERRA.**  The  general  point 
is  P.  Later  on,  P will  be  restricted  to  the  surface. 


Lunar  position  is  specified  in  terms  of  inertial  Cartesian  coordinates,  i = 1, 
2,  3,  by  the  moon  ephemeris  tapes  available  for  use  with  CELEST  and  TERRA. 
Using  one  of  the  coordinate  transformation  computer  routines  available  for  that 
purpose  in  CELEST  and/or  TERRA,  the  inertial  lunar  position  may  be  transformed 
for  any  desired  time  instant  to  the  corresponding  earth-fixed  moon  coordinates  y^^ 


^^Nol  a tide  in  the  strict  sense.  (See  Reference  4,  last  paragraph  of  INTRODUCTION.) 
A new  system  of  computer  programs  for  satellite  geodesy  being  developed  at  NSWC. 


Figure  1.  Position  of  Moon  Relative  to  the  Earth 


Let  now  be  the  moon’s  mass.  G is  Newton’s  constant.  To  find  the  moon’s 
gravitational  potential  about  the  COE,  consider 


r = R - r, 
p L 


Tl  = R-* 


= 


r 

Ml 


GM, 


Ml 


R 


IR- 


(1) 


As  already  stated,  P may  safely  be  assumed  to  be  within  the  close  vicinity  of 
the  earth 


r^  « 50  R 
Ll  > R 
R 6000  km 


(2) 
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Let  7 be  the  angular  distance  between  the  moon  and  point  P.  Let  be  latitude 
and  X be  longitude.  Considering  the  familiar  relationships 


IR  - r^  1 = Ir^  - Rl  = r^  + R^  - 2 rj^  R cos  7 


1 °° 

- 2 — TT  P (cos  7) 

|\-R|  n = or'^^' 

the  potential  of  lunar  mass  attraction  takes  on  the  form 


n-i)  r, 


From  Figure  2 follows 


cos  7 = cos  tJj  cos  i?2  + sin  dj  sin  cos  (Xj  ~ Xj ) 
= sin  sin  + cos  v5,  cos  cos  (X2  " Xj ) 


Also,  it  can  be  shown  that 


P^(cos7)  = P^(sin«p,  )P„(sin^2) 


^ m^=  I (n+“m)!  - ^1 ) (6) 


Finally,  letting  the  lunar  position  correspond  to  index  “1”  and  the  surface  point  to 
index  "2,”  the  expansion  of  the  moon’s  potential  is 


*Wc  intentionally  avoid  distinguishing  between  the  various  types  of  latitude  because  we  expect  that  any  inaccuracy  arising 
from  this  simplification  will  be  masked  by  the  much  larger  uncertainties  which  otherwise  prevail  in  current  tidal  models. 
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Note  that  Equation  (7)  represents  the  moon’s  potential  as  a sum  over  surface 
harmonics. 


V 


n*0 


1 

RH+l 


(8) 


assuming  that  P has  already  been  restricted  to  the  surface  (R  = const),  the  surface 
harmonics  being 


“ a„P  (sin(p)+  2 a Pj"  (sin  tp)  cos  mX 

••  ^ It  —.1  n ni  n 

m - 1 

b„n,P:’(sins^)sinmX 


(9) 


Optionally,  we  may  interpret  Equation  (7)  as  an  expansion  of  the  potential  in  terms 
of  solid  harmonics, 


EI|^(R,(p,X) 


Rn  + I 


(10) 


THE  TIDE  POTENTIAL 

As  it  stands,  0^  is  of  course  not  itself  the  tide  potential  (tide  generating 
potential),  because  the  earth-fixed-reference  system  y^O  of  Figure  1 is  not  an  inertial 
system.  Rather,  the  moon’s  mass  attraction  imparts  to  y*0  the  “guidance” 
acceleration 


Consequently,  any  mass  when  referred  to  y^9 
inertia. 


rL-  (11) 

experiences  the  acceleration,  caused  by 


u - iit 
‘‘in  .2  ‘‘l  • 


Tliis,  in  turn,  can  be  represented  as  the  gradient  of  a potential, 


(12) 
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'•mm 


AIl  ml  a 
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Ml 

a,^  = grad  0,n  = - — (gr^^l  O • 
1*1 


Ml  , _ _ Ml  - 
;^  ‘ • -^L  - - 7 ^ • 

•^L  “^L 


(13) 


For  the  present  purpose,  <^>,1,^  is  of  interest  only  on  the  surface,  r = R,  and 


thus 


0|N  = - — R COS7  - 


- ^ RP,(cos7) 

a 


(14) 


Tlie  total  contribution  of  lunar  mass  attraction  to  the  mechanism  which  causes  the 
tidal  bulge  is  then  the  tide  potential. 


■^L  = '^IN  +'^L 


(15) 


Not  all  of  the  terms  contained  in  Equation  (7)  for  0^^  are  needed  however. 
Closely  following  the  reasoning  of  Reference  8,  Page  70,  0^  was  truncated  after  the 
fourth  term,  resulting  in 


<^L  ^ <'’lO  +<^L1 

Po(cos7)  ^ ^ 
r. 


«L0  = Ml 
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Li  = 


(16) 


R^ 


^L2  = Ml  -J  P2(cos7) 


R^ 

•^LS  = Ml  -T 
*L 
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Thus, 


+'^L0  +^L1  +'^L2  +^3  (17) 

Closer  scrutiny  of  Equation  (17)  shows  that  the  first  cerm  cancels  the  third. 
Also,  the  second  term  is  of  no  consequence  because  its  gradient  vanishes.  The  final 
version  of  the  tide  potential  is  now  quite  simple 

'/'l  3 7)  + ~4  P3(cOS7)  (18) 

1 •'l 

The  Legendre  polynomials  in  this  equation  might  be  specified  from  Equation  (6)  as 
functions  of  the  angular  coordinates  of  both  the  moon  and  the  surface  point. 
However,  because  it  will  greatly  simplify  the  algorithm,  it  is  suggested  to  use 
Equation  (5)  first  to  calculate 


cos  7 = sin  sin  + cos  tp  cos  cos  (X  - ) 


(19) 


and  to  subsequently  evaluate  the  two  Legendre  polynomials. 


TIDE  POTENTIAL  AND  SURFACE  DEFORMATION 

References  7,  8,  and  9 relate  how  elasticity  theory  may  be  applied  to  the 
geophysics  of  the  solid  earth  tide.  Essentially,  the  result  is  as  follows.  Assume  that 
is  one  of  the  harmonic  constituents  of  the  tide  potential  The  radial 

displacement  of  the  earth’s  surface,  caused  by  the  presence  of  this  term,  is 
proportional  to  the  term  itself 


H„(v7,X)  = h^ 


U„(r  = R,X,vJ) 


g 


(20) 


I 


j 
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In  this  equation  occurs  a proportionality  factor  This  is  the  above  mentioned 

geometrical  Love  number  or,  simply,  “Love  number  h^.”  We  regard  it  a constant 
(not  dependent  on  >p  or  X).  Surface  gravity  is  g.  It  is  also  assumed  constant 
(g  go  9.81  m/sec^).  For  further  details,  see  Reference  7,  pp  272  and  273; 
Reference  8,  pp  70-72;  Reference  9,  Chapter  2. 
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Tlie  “tidal  bulge”  may  be  visualized  as  a surface  which  is  displaced  by  a 
distance  H(i;),X)  from  the  unperturbed  earth  surface,  as  illustrated  in  Figure  3.  The 
displacement  H is  the  sum  of  its  components, 


SURFACE  OF  TIDAL  BULGE 


UNPERTURBED  EARTH  SURFACE 


Figure  3.  Geometry  of  Tidal  Bulge 


H(^,X)  = S H^(i^,X). 


The  surface  displacement  now  follows  from  Equations  (21),  (20),  (10),  (8),  (17), 
and  (18). 

u u 

H(^,X)  = hj  — -j  P2(cos7j^)  + h3  — — P3(cos7l)  (22) 

^ S ''l 

where  the  Legendre  polynomials  P,  and  P,  are 


Pjlx)  = -(3x2  _ 


10 


(23b) 


PjCx)  = ^(5x3  - 3x) 
and 

cos7j^  = sin  ip  sin  + cos(p  cosipj^  cos  (X  - X^)  (19) 

POSITION  OF  MOON  AND  TIDE  LAG 

The  three  preceding  sections  demonstrated  how  the  tidal  surface  displacement 
depends,  via  the  tide  potential,  on  the  moon  ephemerides.  The  discussion  was 
entirely  based  on  the  equilibrium  tide  concept  introduced  in  the  LIMITING 
ASSUMPTIONS.  While  thus  in  our  model  the  tidal  elevation  H is  in  phase  with  the 
prevailing  tide  potential  1//^,  there  exists  in  the  real  world  a time  lag  between  H 
and  This  time  lag  has  a twofold  origin.  Firstly,  due  to  the  combined  effects  of 
the  inertia  and  viscosity  inherent  in  the  earth’s  masses,  the  tidal  bulge  will  be 
delayed  with  respect  to  the  moon’s  motion.  At  any  given  time  instant,  the  shape  of 
the  bulge  will  correspond  not  to  the  simultaneous  but  to  a slightly  earlier  position 
of  the  moon  overhead.  And,  secondly,  also  because  of  inertia  and  viscosity,  the 
earth’s  rotation  will  “sweep  along’’  the  bulge,  out  from  under  the  moon,  in  an 
easterly  direction. 

It  would  have  been  quite  impractical  to  account  for  the  complicated 
discrepancy  between  equilibrium  tide  and  actual  tide  in  the  derivations  leading  to  H. 
Instead,  it  was  found  convenient  (Reference  6,  Pages  2 and  9)  to  resolve  it  by 
modifying  the  lunar  position  as  shown  in  Reference  6,  Page  10,  Figure  1.  Thus,  the 
tidal  retardation  may  be  reflected  in  our  model  by  introducing  a time  delay  into 
the  lunar  ephemeris.  The  rotational  advance  of  the  bulge  may  be  introduced  by 
suitably  rotating  the  orbital  plane  associated  with  the  motion,  relative  to  the  earth, 
of  the  moon. 

No  further  mathematical  details  appear  necessary.  The  equations  which  reflect 
the  tidal  lag  in  the  second  paragraph  of  the  next  section  on  the  computer  algorithm 
will  illustrate  the  matter  sufficiently.  But  it  still  ought  to  be  mentioned  that  there 
is  a mathematical  flaw  in  this  part  of  the  algorithm.  Namely,  the  algorithm  first 
retards  the  moon  by  setting  back  the  inertial  lunar  coordinates  to  a sliglitly  earlier 
orbital  position.  Subsequently  it  computes  the  earth-fixed  coordinates  of  the  moon. 
The  advance  of  the  lunar  orbit  plane  is  then  accomplished  by  suitably  manipulating 
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the  earth-fixed  moon  coordinates.  The  author  is  obliged  to  Drs.  C.  Oesterwinter  and 
J.  OToole  of  the  Astronautics  and  Geodesy  Division  for  pointing  out  to  him  that 
this  procedure  is  not  invariant  with  respect  to  the  epoch  of  the  lunar  ephemeris  as 
any  procedure  of  this  type  should  be. 

Let  X‘  be  the  lunar  position  vector  at  time  instant  t,  referred  to  the  inertial 
frame  corresponding  to  time  t.  Let  Y'  be  the  corresponding  earth-fixed  position  of 
the  moon.  Let  the  matrix  product*  ABCD  be  the  transformation  matrix  leading 
from  the  inertial  frame  to  the  earth-fixed  frame, 


Y'  = ABCD  X‘  . (24) 

D is  the  Precession  Matrix,  leading  from  the  mean  equator  and  equinox  of  the 
epoch  frame  to  that  of  date.  The  Nutation  Matrix  is  C which  continues  the  just 
mentioned  rotation  on  to  the  frame  asssociated  with  the  true  equator  and  equinox 
of  date.  B continues  the  transformation  to  the  earth-fixed  frame  of  rotation,  and  A 
transforms  to  the  frame  of  figure  (correction  for  polar  motion). 

Now,  let  “40”  and  “50”  indicate  two  different  epochs,  say  1940.0  and  1950.0, 
each  defining  a particular  inertial  frame.  Then, 


X'  = ,X'» 

,„X*» 

X'  - Dso,,D40,soX'" 


(25) 


and,  defining  L to  be  a matrix  representing  the  position  difference  between  the 
ephemeris  moon  and  the  fictitious  moon  (or  part  of  that  position  differential,  such 
as  the  retardation  of  the  moon  in  its  orbit). 


Y'  = ABCD5„_,(LX^«)  = ABCDjo^.LD.o^j^X^" 


ABCD,o^,LX^» 


(26) 


« 

From  unpublished  working  notes  by  W.  Groeger,  on  the  Development  of  the  Equations  of  Motion  for  the  TERRA 
Computer  Program  System,  August  1973. 
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as,  generally,  L and  D, 


will  not  commute 


0,50 


^ ^4  0,5  0 *^4  0,5  0*- 


(27) 


However,  everyone  agrees  with  the  author  that  the  tide  lag  algorithm  appearing 
below  in  the  COMPUTER  ALGORITHM  section,  because  of  the  small  magnitude  of 
the  effect,  is  a legitimate  approximation  to  reality  despite  its  mathematical 
inconsistency. 


CONTRIBUTION  BY  THE  SUN 

Everything  said  in  the  above  sections  may  be  equally  applied  to  the 
deformation  of  the  earth  by  the  tidal  action  of  the  sun.  Note  especially  that  the 
Love  numbers  h^  are  geophysical  constants  and  depend  on  the  structure  of  the 
earth  alone.  They  consequently  apply  equally  to  tlie  lunar  tide  and  to  the  solar 
tide. 


The  sun  is  much  more  distant  from  the  earth  than  the  moon. 


(28) 


suggesting  that  only  the  first  term  be  considered  in  the  equation  for  the  surface 
displacement. 


W 2 

H((p,X)  = hj  — -jP2(cos75) 

8 Jc 


(29) 


cos  75  = sin  tp  sin  + cos  yp  cos  tp^  cos  (X  - X^ ) 


(30) 
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COMPUTER  ALGORITHM 


I 


L 


INPUT  DATA 

Geophysical  Constants 

g Approximation  of  earth’s  surface  gravity,  in  terms  of  m/sec^.  Suggested 

value:  g = 9.81  m/sec^. 

R Semimajor  axis,  in  meters,  of  a suitable  reference  ellipsoid.  Use  a value 

compatible  with  the  data  for  the  particular  CELEST  run.  Or  else, 
R = 6 378  150  m ought  to  be  satisfactory. 

hj  Love  number.  Suggested  value  for  trial  runs:  hj  = 0.6. 

Moon’s  gravitational  constant,  in  cm^/sec^.  = G where  Mj^  is  the 
lunar  mass.  Use  a value  compatible  with  CELEST  input.  Suggested  value  for 
trial  runs:,  = 4.9177*  10*®  cm^/sec^. 

Pg  Sun’s  gravitational  constant,  in  terms  of  cm^/sec^.  Pg  = G Mg  where  Mg  is 

the  solar  mass.  Use  a value  compatible  with  the  CELEST  input.  Suggested 
value  for  trial  runs:  Pg  = 1.3291  * 10^*  cm®/sec^. 

Station  Coordinates 

<fi  Latitude  of  station  (geocentric  or  geodetic  etc.).  May  be  expected  to  be 
available  in  terms  of  degrees,  from  a station  catalog.  Express  this  either  in 
decimal  fraction  of  degrees  or  in  radians,  to  be  compatible  with  those 
trigonometric  functions  where  <p  is  the  argument. 

X Station  longitude.  The  comments  just  made  for  the  station  latitude  apply 

here  also. 


Position  of  Moon  and  Sun 


t Time  instant  for  which  the  tidal  surface  displacement  H is  to  be  found. 
Specify  t in  terms  of  year,  day  UT,  and  seconds  elapsed  since  beginning  of 
day  UT. 


j 


At  Tidal  lag  parameter  in  seconds. 

tj  Earth’s  sidereal  rate  of  rotation, 

u = 4.178  074  622  ♦10'^  deg/sec 
= 7.  292  115  855  ♦ lO'^  rad/sec 

Have  available  suitable  ephemerides  for  moon  and  sun,  such  as  the  TERRA/CELEST 
moon  and  sun  ephemeris  tapes,  plus  the  associated  interpolation  routine. 

When  to  Evaluate  the  Tidal  Displacement  During  the  Data 
Blocks  on  the  CELEST  Observation  Tape 

tp^  Time,  in  seconds,  of  closest  approach  (of  satellite  to  station)  during  data  block 
(satellite  pass)  on  CELEST  observation  tape. 

tj  First  observing  time  value,  in  seconds,  for  the  individual  data  block  on  the 

CELEST  observing  tape. 

T Time  increment,  in  seconds,  necessary  to  specify  the  time  instants  for  which 

the  tidal  displacement  is  to  be  found. 

COLLECTION  OF  EQUATIONS  FOR  CODING 
Position  of  the  Moon 

For  the  time  instant  (t  - At),  obtain  the  inertial  Cartesian  coordinates  of  the 

moon,  x^,  by  interpolation  from  the  moon  ephemeris  tape.  Note  that  a 

suitable  interpolation  routine  exists  for  use  with  the  tape. 

Transform  the  latter  coordinates  to  the  corresponding  earth-fixed  Cartesian 
coordinates,  x'^,  y'^,  z'^.  The  necessary  transformation  procedure  is  readily  available 
in  CELEST  and/or  TERRA  as  a subroutine.  Reference  is  also  made  to  the 
POSITION  OF  MOON  AND  TIDE  LAG  section.  If  not  already  expressed  in  terms 
of  kilometers,  now  write  x'^,  y'^  and  in  terms  of  this  unit. 
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Figure  4.  Schematic  of  the  Algorithm 


16 


Now  transform  to  earth-fixed  spherical  coordinates. 


! 

5 

u 


Zi. 

= tan' ' ^ 

+V(x'l)2  +(y'j^)2 

X'l  = tan'’ 

L, 


(lOlA) 

(lOlB) 

(lOlC) 


taking  the  customary  precautions  to  assure  that  the  resulting  angle  values  are  located 
in  the  correct  quadrant. 

Then  advanc.'  the  position  of  the  moon  in  terrestrial  longitude  (longitude  east) 
to  obtain  the  earth-fixed  position  of  the  fictitious  (tide  generating)  moon 


(102  A) 

= < 

(102B) 

= Xl  +wAt 

(I02C) 

Sun  Position 

Now  calculate  the  position  of  the  fictitious  sun  (tide  generating  sun).  Code  for 

this  purpose  a procedure  similar  to  that  just  specified  for  the  moon. 

Station  Coordinates 

The  station  coordinates  can  be  expected  to  be  read  off  a catalog  for  the 

observing  stations  directly  in  terms  of  longitude  and  latitude.  For  the  present 
purpose,  geodetic,  geographic,  and  geocentric  latitude  may  be  freely  interchanged  and 
will  in  the  following  be  referred  to  as  “latitude.”  It  may  reasonably  be  assumed 

that  any  arbitrary  surface  points,  other  than  observing  stations,  will  also  be  specified 
in  terms  of  their  latitudes  and  longitudes. 


liL 
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Computation  of  Surface  Displacement 

Now  find  the  surface  displacement,  H,  associated  with  a particular  station  or 
surface  point  as  well  as  with  a specific  time  instant  t.  The  radial  displacement  is  H. 

A positive  value  of  H implies  that  the  station  has  increased  its  distance  from  the 

j COE.  Tidal  movement  of  the  station,  parallel  to  the  surface,  has  been  neglected  for 

I the  purpose  of  the  analysis  which  led  to  the  present  algorithm.  From  the  equations 

below,  H wDl  result  in  terms  of  meters. 

I 
1 

u.  ii 

H((p,X,t)  = 10"'^hj—  3 P,(cos7,  )+ lO^'^hj  — 3 P-fcos-yg) 

2 g 2 L 2 g ^ 2 s 

Pjfx)  = ^(3x2-  I) 

cos  'P  *>*1  "Pl  (105) 

cos  75  = sin  \p  sin  <^5  + cos  cos  ipg  cos  (X  - Xj ) (106) 

When  to  Evaluate  the  Tidal  Displacement  During  the  Data 
Blocks  on  the  CELEST  Observation  Tape 

First  Option  For  each  data  block  on  the  observation  tape  (satellite  pass), 

evaluate  H at  the  time  of  closest  approach  (of  the  satellite  to 
the  station)  H = H(tp^).  Regard  this  as  the  value  of  H 
throughout  the  satellite  pass. 

Second  Option  For  each  satellite  pass  (data  block  on  the  observation  tape), 

evaluate  H for  the  first  observing  time  value,  t, , and, 
subsequently  for  tj  = t,  +t,  tj  = tj  +t,  ...»  t^^^j  = t^^  +t,  ...  , 
up  to  the  end  of  the  satellite  pass. 

For  the  observing  time  t,  t^^  < t < t^^i,  assign  to  H the  value 
H(t„). 

T is  supposed  to  be  specified  by  ‘the  CELEST  user  via  the 
program  input.  Make  provision  for  the  user  to  insert  the  chosen 
value  of  r together  with  the  input  data  for  the  individual 
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L- A 


(103) 

(104) 


computer  run.  If  the  present  second  option  is  chosen  at  all,  it  is 
quite  likely  that  the  value  r = 10  min  (600  sec)  will  be  selected 
because  this  appears  to  be  a suitable  value  for  a wide  range  of 
satellite  altitudes  (pass  lengths  between  20  min  and  10  hr). 


COMPUTER  PROGRAM  OUTPUT 

It  can  be  expected  that  the  present  computer  program  will  serve  as  a 
subroutine  during  CELEST  runs.  Whenever  required,  CELEST  will  specify  latitude 
and  longitude  for  a surface  point  plus  a time  value.  The  present  routine  will  then 
calculate  the  surface  displacement  associated  with  the  station  position  and  the  time 
instant. 

Optionally,  provision  may  be  made  to  input  for  any  given  station  position,  a 
table  of  time  values  and  to  automatically  find  all  associated  surface  displacements. 
The  latter  procedure  is  expected  to  be  useful  for  the  program  checkout  and  also  for 
certain  geophysical  studies  made  independently  of  CELEST. 


QUANTITATIVE  RESULTS  FOR  THREE  TYPICAL  STATION  LOCATIONS 


The  numerical  checkout  of  the  just  specified  computer  program  was  based  on 
the  input  parameters  listed  in  the  first  paragraph  of  the  COMPUTER  ALGORITHM. 
Surface  displacements  (total  value  H in  centimeters,  lunar  component  HMOON  in 
centimeters  and  solar  component  HSUN  in  centimeters)  were  computed  for  each  of 
three  typical  stations  at  ten-minute  intervals  for  an  entire  day  (Day  UT  No.  88  of 
Year  1977). 
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Table  1.  Surface  Displacement  and  Lunar  and  Solar  Components  of 
Surface  Displacement  for  Station  No.  1 

(Station  longitude  = 0 deg,  station  latitude  = 0 deg,  HMOON,  HSUN,  and  H are  in 
centimeters.  Time  is  listed  in  terms  of  seconds  from  Midniglit  UT  of  Day  88  of 
Year  1977.) 


TIME 

■= 

a 

SEC 

HrtOON 

-8.05 

HSUN 

9.86 

H 

£ 

1.82 

TIME 

600 

SEC 

HMOON 

-8  .59 

HSUN 

= 

9.87 

H 

£ 

1.28 

TIME 

= 

1200 

SEC 

HMCCN 

-9.04 

ISUN 

= 

9.82 

H 

£ 

.78 

TIME 

IflCJ 

SEC 

HMCCN 

-9.41 

HSUN 

9.72 

H 

£ 

.32 

TIME 

= 

2400 

SEC 

HMOON 

= 

-9.6  f 

HSUN 

9.56 

H 

£ 

-.11 

TIME 

30  00 

SEC 

HMCCN 

-9.85 

HSUN 

9.35 

H 

= 

-.50 

TIME 

36  00 

SEC 

HMOON 

-9.93 

HSUN 

= 

9.08 

H 

= 

-.84 

TIME 

4200 

SEC 

HMOCN 

-9.93 

HSUN 

8.77 

H 

-1.1-4 

TIME 

48  30 

SEC 

HMOON 

= 

-9.79 

HSUN 

= 

8.40 

H 

£ 

-1.39 

TIME 

54  00 

SEC 

HMCON 

-9.5  < 

HSUN 

8.00 

H 

£ 

-1.58 

TIME 

s 

6000 

SEC 

HMCCN 

-9.27 

fSUN 

■= 

7.54 

H 

= 

-1.73 

TIME 

= 

6600 

SEC 

HMCCN 

= 

-8.87 

HSUN 

= 

7.05 

H 

£ 

-1.81 

time 

72C0 

SEC 

HMOON 

-8.21 

HSUN 

6.5? 

H 

£ 

-1.  85 

TIME 

= 

7800 

SEC 

HMCCN 

-7.  er 

HSUN 

= 

5.97 

H 

£ 

-1.83 

TIME 

8400 

SEC 

HMOON 

-7.14 

HSUN 

= 

5.39 

H 

£ 

-1.  75 

time 

9000 

SEC 

HMCON 

= 

-6.4( 

HSUN 

= 

4.78 

H 

£ 

-1.61 

time 

9600 

SEC 

HMCCN 

= 

-5  .5  1 

HSUN 

= 

4.16 

H 

£ 

-1.42 

time 

= 

10200 

SEC 

HMCON 

= 

-4.71 

HSUN 

= 

3.53 

H 

£ 

-1.18 

TIME 

r 

10800 

SEC 

HMCON 

-3.76 

HSUN 

2.88 

H 

£ 

-.88 

time 

r 

11400 

SEC 

HMCCN 

= 

-2.77 

HSUN 

= 

2 .24 

H 

£ 

-.53 

TIME 

= 

120C0 

SEC 

HMCCN 

-1.73 

HSUN 

s 

1.59 

H 

£ 

-.14 

TIME 

= 

12600 

SEC 

HMCCN 

- . 64 

fSUN 

= 

.95 

H 

£ 

.31 

time 

= 

1 3200 

SEC 

HMCCN 

= 

.47 

HSUN 

= 

.32 

H 

£ 

.80 

TIME 

= 

1 38  00 

SEC 

HMCCN 

= 

1.61 

HSUN 

£ 

-.29 

H 

£ 

1.32 

TIME 

= 

14403 

SEC 

HMCCN 

2.77 

HSUN 

= 

-.88 

H 

£ 

1.89 

TIME 

15000 

SEC 

HMCCN 

= 

3.93 

tSUN 

-1.44 

H 

£ 

2.49 

TIME 

156C0 

SEC 

HMOON 

■= 

5 .10 

HSUN 

£ 

-1.98 

H 

£ 

3.12 

TIME 

= 

16200 

SEC 

HMCCN 

6.25 

HSUN 

= 

-2.48 

H 

£ 

3.77 

TIME 

s 

16600 

SEC 

HMCCN 

= 

7,3* 

HSUN 

£ 

-2.94 

H 

£ 

4.45 

TIME 

= 

174C0 

SEC 

HMCCN 

= 

8.5  3 

HSUN 

£ 

-3.37 

H 

£ 

5.14 

time 

r 

1 8C00 

SEC 

HMCCN 

r 

9.5* 

HSUN 

£ 

-3.74 

H 

£ 

5.85 

TIME 

= 

18600 

SEC 

HMCCN 

■= 

10.63 

HSUN 

£ 

-4.08 

H 

£ 

6.56 

TIME 

z 

192C0 

SEC 

HMOON 

r 

11.6? 

tSUN 

£ 

-4.36 

H 

£ 

7.27 

TIME 

= 

19800 

SEC 

HilOCN 

12.57 

HSUN 

£ 

-4.59 

H 

£ 

7.98 

TIME 

s 

2C4C0 

SEC 

HMCON 

= 

13.45 

HSUN 

s 

-4.76 

H 

£ 

8.68 

TIME 

r 

210C0 

SEC 

HMCON 

= 

14. 2t 

HSUN 

£ 

-4.88 

H 

£ 

9.37 

TIME 

= 

21600 

SEC 

HMCCN 

= 

14  .99 

fSUN 

£ 

-4.95 

H 

£ 

10.05 

TIME 

s 

22200 

SEC 

HMCON 

= 

15.65 

HSUN 

£ 

-4.96 

H 

£ 

10.69 

TIME 

= 

228  00 

SEC 

HMOON 

16.23 

HSUN 

£ 

-4.91 

H 

= 

11.32 

t 

1 
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Table  1.  Surface  Displacement  and  Lunar  and  Solar  Components  of 
Surface  Displacement  for  Station  No.  1 (Continued) 


TIME 

23i»ca 

SEC 

HMCCN 

= 

16.71 

HSUN 

X 

-4.81 

H 

X 

11.91 

TIME 

= 

2<«oaa 

SEC 

HMCON 

X 

17.11 

»SUN 

X 

-4.65 

H 

X 

12.46 

TIME 

s 

24600 

SEC 

HMCON 

X 

17.41 

HSO'N 

X 

-4  .43 

H 

X 

12.98 

TIME 

= 

25200 

SEC 

HMCON 

X 

17.61 

HSUN 

X 

-4.17 

H 

X 

13.45 

TIME 

s 

25800 

SEC 

HMCCN 

X 

17.72 

HSUN 

X 

-3.85 

H 

X 

13.87 

TIME 

s 

26400 

SEC 

HMCON 

X 

17.73 

6SUN 

X 

-3.49 

H 

X 

14.24 

time 

s 

27000 

SEC 

HMCCN 

X 

17. 6t 

HSUN 

X 

-3.08 

H 

= 

14.5b 

TIME 

= 

2 7600 

SEC 

HMCCN 

X 

17.95 

HSUN 

X 

-2.63 

H 

X 

14.82 

TIME 

= 

28200 

SEC 

HMCCN 

X 

17.11 

HSUN 

X 

-2.14 

H 

X 

15.  03 

TIME 

= 

28800 

SEC 

HMCCN 

X 

16.79 

► SUN 

X 

-1.61 

H 

X 

15.17 

time 

29400 

SEC 

HnOCN 

X 

16.3] 

HSUN 

X 

-1.06 

H 

X 

15.26 

TIME 

= 

30  C 00 

SEC 

HMCCN 

X 

15.75 

HSUN 

X 

-.47 

h 

X 

15.28 

TIME 

30600 

SEC 

HMCCN 

X 

15.10 

HSUN 

X 

.13 

H 

X 

15.23 

TIME 

r 

31200 

SEC 

HMCON 

X 

14.38 

(SUN 

X 

.75 

H 

X 

15.13 

TIME 

= 

31800 

SEC 

HMCON 

X 

13.57 

HSUN 

X 

1.39 

H 

X 

14.96 

TIME 

32400 

SEC 

HMCON 

X 

12.70 

HSUN 

X 

2.03 

H 

X 

14.73 

TIME 

33000 

SEC 

HMCON 

X 

11.77 

HSUN 

X 

2.68 

H 

X 

14.45 

TIME 

= 

33600 

SEC 

HMCON 

X 

10.77 

HSUN 

X 

3.32 

H 

X 

14.10 

TIME 

= 

34200 

SEC 

HMCON 

X 

9.73 

HSUN 

X 

3.96 

H 

X 

13.69 

TIME 

= 

34800 

SEC 

HMCCN 

X 

8.65 

HSUN 

X 

4.59 

H 

X 

13.24 

TIME 

= 

35400 

SEC 

HMCCN 

X 

7.53 

HSUN 

X 

5.2C 

H 

X 

12.  73 

TIME 

= 

36003 

SEC 

HMCON 

X 

6.36 

HSUN 

X 

5.79 

H 

X 

12.17 

TIME 

= 

36600 

SEC 

HMCON 

X 

5.22 

HSUN 

X 

6.35 

H 

X 

11.57 

TIME 

= 

37200 

SEC 

HMCCN 

X 

4.04 

HSUN 

X 

6.89 

H 

X 

10.93 

TIME 

= 

37800 

SEC 

HMCON 

X 

2 .86 

HSUN 

X 

7.39 

H 

X 

10.25 

TIME 

38403 

SEC 

HMCON 

X 

1.69 

HSUN 

X 

7.85 

H 

X 

9.55 

TIME 

= 

39000 

SEC 

HMCON 

X 

.54 

HSUN 

X 

8.28 

H 

X 

8.81 

TIME 

396C0 

SEC 

HMCCN 

X 

- .60 

HSUN 

X 

8.65 

H 

X 

6.06 

TIME 

= 

40200 

SEC 

HMCCN 

X 

-1.70 

HSUN 

X 

8.98 

H 

X 

7.29 

TIME 

= 

408CO 

SEC 

HMCCN 

X 

-2.76 

HSUN 

X 

9.26 

H 

X 

6.50 

TIME 

= 

41400 

SEC 

HMCON 

X 

-3.77 

HSUN 

X 

9.49 

H 

X 

5.72 

TIME 

s 

42000 

SEC 

HMCCN 

X 

-4  .74 

HSUN 

X 

9.b7 

H 

X 

4.93 

TIME 

= 

42600 

SEC 

HMCON 

X 

-5.64 

HSUN 

X 

9.79 

H 

X 

4.15 

TIME 

r 

43200 

SEC 

HMCCN 

X 

-6.47 

HSUN 

X 

9.85 

H 

X 

3.38 

TIME 

= 

43800 

SEC 

HMCON 

X 

-7.23 

HSUN 

X 

9.86 

H 

X 

2.63 

time 

= 

44400 

SEC 

HMCCN 

X 

-7.91 

HSUN 

X 

9.81 

H 

X 

1.90 

TIME 

s 

45C0G 

SEC 

HMCCN 

X 

-8.51 

HSUN 

r 

9.71 

H 

X 

1.20 

TIME 

= 

456C0 

SEC 

HMCCN 

X 

-9. Cl 

HSUN 

X 

9.55 

H 

X 

.53 

TIME 

= 

46200 

SEC 

HMCON 

X 

-9.43 

HSUN 

X 

9.34 

H 

X 

-.10 

TIME 

s 

46800 

SEC 

HMCCN 

X 

-9.75 

HSUN 

X 

9.07 

H 

X 

-.68 

TIME 

= 

47400 

SEC 

HMCCN 

X 

-9.98 

HSUN 

X 

8.75 

H 

X 

-1.22 

TIME 

= 

48000 

SEC 

HMCCN 

X 

-10. 1C 

HSUN 

X 

8.39 

H 

X 

-1.71 

TIME 

= 

48600 

SEC 

HMCCN 

X 

-10.13 

HSUN 

X 

7.98 

h 

X 

-2.15 

TIME 

s 

49230 

SEC 

HMCCN 

X 

-10.05 

HSUN 

X 

7.53 

H 

X 

-2.52 

TIME 

= 

49800 

SEC 

HMCCN 

X 

-9.ee 

HSUN 

X 

7.04 

H 

X 

-2.84 

TIME 

= 

50400 

EEC 

HMCON 

X 

-9.6( 

HSUN 

X 

6.51 

H 

X 

-3.09 

TIME 

X 

51000 

SEC 

HMCCN 

X 

-9.23 

HSUN 

X 

5.96 

H 

X 

-3.27 

TIME 

= 

5 1600 

SEC 

HMCON 

X 

••8.76 

HSUN 

X 

5.37 

H 

X 

-3.39 

TIME 

X 

52200 

SEC 

HMCCN 

X 

-8.21 

HSUN 

X 

4.77 

H 

X 

-3.44 

TIME 

s 

52800 

SEC 

HMCCN 

X 

-7.56 

HSUN 

X 

4.15 

H 

X 

-3.42 

TIME 

X 

53400 

SEC 

HMCON 

X 

-6.84 

HSUN 

X 

3.51 

H 

X 

-3.32 

TIME 

X 

54000 

SEC 

HMCCN 

X 

-6.  C3 

HSUN 

X 

2.87 

H 

X 

-3.16 

TIME 

= 

54600 

SEC 

HMCON 

X 

-5.15 

HSUN 

X 

2.22 

H 

X 

-2.93 

TIME 

X 

55200 

SEC 

HMCON 

X 

-4.21 

HSUN 

X 

1.58 

H 

X 

-2.63 
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Table  1.  Surface  Displacement  and  Lunar  and  Solar  Components  of 
Surface  Displacement  for  Station  No.  1 (Continued) 


TIME 

s 

55800 

SEC 

HMCOM 

2 

-3.31 

HSUN 

2 

.9^ 

h 

2 

-3.37 

TIME 

= 

561.00 

SEC 

H'1CC9 

2 

- 3 . 1 E 

HSUN 

2 

. 31 

h 

2 

-1.84 

time 

57  0 CO 

SEC 

HMCC.'I 

2 

- 1 . 0 E 

HSUN 

2 

-.3C 

H 

Z 

-1.35 

TIME 

= 

5 76  0 0 

SEC 

HdCCN 

2 

.CS 

HSUN 

2 

- .89 

H 

Z 

-.3  0 

TIME 

= 

5 e3C0 

SEC 

HdCC.N 

2 

1.35 

HSUN 

2 

-1.45 

H 

= 

-.31 

TIME 

= 

5 8800 

SEC 

HMOCN 

2 

2 .44 

^ SUN 

2 

-1.99 

h 

2 

.45 

TIME 

501.00 

SEC 

HMOCd 

2 

3.64 

HSU  iN 

2 

-3.49 

h 

2 

1.15 

TIME 

= 

b 0 0 : 0 

SEC 

HMCCN 

2 

4.84 

HSUN 

2 

-3.95 

h 

2 

1.8  9 

TIME 

= 

6 0600 

SEC 

HMCCN 

2 

6.04 

HSUN 

2 

-3.37 

H 

= 

3.87 

TIME 

61200 

SEC 

HMCCN 

2 

7.23 

1 SUN 

2 

-3.7j 

ti 

2 

3. -48 

TIME 

= 

6 18  00 

SEC 

HMCOM 

2 

8.39 

)’3UN 

= 

-4.0  8 

H 

= 

4.31 

TIME 

s 

62i.CC 

SEC 

HMOCN 

2 

9 .53 

HSUN 

2 

-4.36 

h 

2 

5.16 

TIME 

= 

6 3 0 00 

SEC 

HMCCN 

2 

10.83 

HSUN 

2 

-4  .59 

H 

2 

6.03 

TIME 

= 

636  00 

SEC 

HMCCN 

2 

11.67 

)SUN 

2 

-4 . 7t 

H 

2 

6.91 

TIME 

61.200 

SEC 

HMCCN 

2 

13.68 

H SUN 

2 

-u  .88 

H 

2 

7.78 

TIME 

= 

648  CO 

SEC 

HMCCN 

2 

1 3 .bC 

HSUN 

2 

-4.95 

h 

2 

8.63 

TIME 

654  CO 

SEC 

HMCCN 

2 

14.47 

hSCN 

2 

-4.95 

fi 

r 

9.53 

TIME 

= 

6 6 0 jC 

SEC 

HMCC.^ 

2 

15.37 

HSUN 

2 

-4.9  1 

H 

2 

i:.36 

TIME 

66600 

SEC 

HMOCN 

2 

15.98 

HSUN 

2 

-H.ao 

H 

2 

11.18 

TIME 

= 

672C0 

SEC 

MM  CCN 

2 

16  .63 

HSUN 

2 

-4.84 

H 

2 

11.98 

TIME 

= 

67PCQ 

SEC 

HMCCN 

2 

17.18 

)SUN 

2 

-4.43 

h 

2 

13.7  3 

TIME 

684C0 

SEC 

HMCON 

2 

17.6  : 

HSUN 

2 

-4. IE 

H 

= 

1 3 . s 5 

TIME 

6SCu0 

SEC 

H -ICON 

2 

1 7 .97 

HSUN 

2 

- 3.83 

h 

2 

14.12 

TIME 

= 

69600 

SEC 

HMCCN 

2 

18.33 

H5LN 

2 

-3.48 

M 

2 

1m  . 75 

TIME 

r 

7 0200 

SEC 

HM  CCN 

2 

18  .38 

ISUN 

2 

-3. 07 

h 

2 

15.31 

TIME 

= 

7 0 8 0 0 

SEC 

HMCCN 

2 

18.44 

H SUN 

= 

-3  .63 

H 

2 

15.83 

TIME 

z 

71400 

SEC 

H M C G N 

2 

18  .3  ‘ 

HSUN 

2 

-3.13 

2 

16. 3o 

TIME 

= 

720  CO 

SEC 

HMCCN 

2 

18.34 

hSLN 

= 

-1.61 

H 

2 

16.6m 

TIME 

= 

7 26  C0 

SEC 

HMCCN 

2 

1 7 .9E 

8 SUN 

= 

-1 .05 

h 

2 

16.95 

TIME 

= 

7 3200 

SEC 

HMCON 

2 

1 7 .6E 

HSUN 

2 

-.47 

h 

2 

17.18 

TIME 

2 

7 38  00 

SEC 

HMCCN 

2 

17.30 

HSUN 

2 

.14 

H 

2 

17.34 

TIME 

X 

744C0 

sec 

HMCON 

2 

16.66 

HSC'N 

2 

. 76 

H 

2 

17.42 

TIME 

= 

7 EC  00 

SEC 

HMCCN 

2 

16.03 

)SCN 

2 

1.40 

H 

2 

17.43 

TIME 

2 

756  0 J 

SEC 

HM  CG  . 

2 

15.3.' 

H SUN 

2 

3 .C4 

H 

2 

17.36 

TIME 

= 

76200 

SEC 

HMCON 

r 

14.53 

HSUN 

2 

3.68 

h 

2 

17.31 

TIME 

2 

76800 

SEC 

HMCCN 

2 

13. 6t 

HSCN 

2 

0.33 

H 

2 

10.99 

TIME 

2 

7 74  C O 

SEC 

HMCCN 

2 

13.73 

HSUN 

2 

3.97 

H 

2 

10.69 

TIME 

2 

7 8 000 

SEC 

HMCON 

2 

11.73 

HSUN 

2 

4.59 

h 

2 

16.31 

TIME 

2 

7 86  00 

SEC 

HMCCN 

2 

1 : .67 

HSUN 

2 

5.30 

M 

2 

15.97 

TIME 

2 

7920J 

SEC 

HMCCN 

2 

9.57 

8 Son 

2 

5.79 

H 

2 

15.36 

TIME 

2 

7 ESCO 

sec 

HMCON 

2 

8 .43 

HSUN 

2 

6.36 

H 

= 

14.78 

TIME 

= 

8 04  00 

SEC 

HMCCN 

2 

7.35 

HSCN 

2 

6.89 

h 

2 

14.1m 

TIME 

2 

81GC0 

SEC 

HMCGN 

2 

6.0) 

HSLN 

2 

7.39 

H 

2 

13.45 

TIME 

2 

8 1600 

SEC 

HMCON 

2 

4.3E 

HSUN 

2 

7.85 

H 

2 

13. 7 j 

TIME 

2 

8 22''J 

SEC 

HM  CGN 

2 

3.63 

HSUN 

2 

8.37 

h 

2 

11.90 

TIME 

2 

82800 

SEC 

HMCGN 

2 

3.43 

HSLN 

2 

8.65 

H 

2 

11.07 

TIME 

2 

8 34  CO 

sec 

HMCON 

2 

1.21 

8 SUN 

2 

8.98 

h 

2 

10.19 

TIME 

2 

8 4 0 0 0 

SEC 

HM  CON 

2 

.0  3 

HSUN 

2 

9.3t 

H 

2 

9.39 

TIME 

2 

848  CO 

SEC 

HMCGN 

2 

-1.13 

HSUN 

2 

9.49 

h 

2 

8.36 

TIME 

2 

85200 

SEC 

HMCO  N 

2 

-3.34 

HSUN 

2 

9.66 

h 

= 

7.43 

TIME 

2 

85  8 CO 

SEC 

HMCCN 

2 

-3.31 

8SUN 

2 

9.78 

H 

2 

6.47 
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Table  2.  Surface  Displacement  and  Lunar  and  Solar  Components  of 
Surface  Displacement  for  Station  No.  2 

(Station  longitude  = 0 deg,  station  latitude  = 30  deg,  HMOON,  HSUN,  and  H are 
in  centimeters.  Time  is  listed  in  terms  of  seconds  from  Midnight  UT  of  Day  88  of 
Year  1977.) 


TIME 

Q 

SEC 

HMOCN 

= 

-6.04 

HSUN 

5.44 

H 

z 

-.60 

TIME 

600 

SEC 

HMCCN 

= 

-6,74 

HSUN 

5.45 

H 

z 

-1.29 

TIME 

= 

12C0 

SEC 

HMCCN 

= 

-7.3t 

HSLN 

= 

5 .41 

H 

z 

-1.97 

TIME 

18C0 

SEC 

HMOON 

-7.95 

fSUN 

s 

5.34 

H 

z 

-2.61 

TIME 

24C0 

SEC 

HMOCN 

-a  .45 

HSUN 

= 

5.22 

H 

Z 

-3.23 

TIME 

= 

3000 

SEC 

HMOON 

= 

-8.88 

HSUN 

= 

5.06 

H 

Z 

-3.81 

TIME 

= 

36C0 

SEC 

HMCCN 

= 

-9.24 

HSUN 

z 

4.87 

H 

z 

-4.36 

TIME 

<«2C0 

SEC 

HMCON 

= 

-9.52 

ISUN 

= 

4.64 

H 

Z 

-4.08 

TIME 

= 

L8CU 

SEC 

HMOON 

= 

-9.72 

HSUN 

z 

4.38 

h 

z 

-5.36 

TIME 

= 

5400 

SEC 

HMOON 

= 

-9.87 

HSUN 

z 

4.08 

H 

z 

-5.79 

TIME 

= 

6000 

SEC 

HMCCN 

= 

-9.94 

HSUN 

Z 

3.7t> 

H 

Z 

-6.18 

TIME 

= 

6600 

SEC 

HMCCN 

= 

-9.93 

l-SUN 

z 

3.4C 

H 

Z 

-6.53 

TIME 

72C0 

SEC 

HMOON 

= 

-9.85 

HSUN 

= 

3.02 

H 

z 

-6.83 

TIME 

= 

7800 

SEC 

HMCON 

= 

-9.71 

HSUN 

Z 

2,62 

H 

= 

-7.09 

TIME 

T 

84C0 

SEC 

HMCON 

= 

-9.5C 

HSUN 

Z 

2.20 

H 

Z 

-7.30 

TIME 

= 

9000 

SEC 

HMCCN 

= 

-9.22 

fSUN 

z 

1.76 

H 

z 

-7.46 

TIME 

= 

9600 

SEC 

HMCCN 

= 

-8.89 

HSUN 

z 

1.31 

H 

Z 

-7.57 

TIME 

= 

1C200 

SEC 

HMOON 

= 

-8  .50 

HSUN 

z 

.86 

H 

= 

— 7.64 

TIME 

= 

1 08  00 

SEC 

HMCCN 

= 

-8. Of 

HSUN 

z 

.40 

H 

z 

-7.66 

TIME 

s 

11400 

SEC 

HMCCN 

= 

-7  .57 

fSUN 

z 

-.07 

H 

= 

-7.63 

TIME 

= 

12C0O 

SEC 

HMCON 

-7  .Ct 

HSUN 

z 

-.53 

H 

Z 

-7.56 

TIME 

12600 

SEC 

HMCON 

= 

-6.47 

HSUN 

z 

-.98 

H 

z 

-7.45 

TIME 

= 

1 3200 

SEC 

HMCCN 

= 

-5.87 

HSUN 

z 

-1.43 

H 

z 

-7.30 

TIME 

= 

1 3800 

SEC 

HMCON 

= 

-5.24 

(SUN 

z 

-1.86 

H 

= 

-7.10 

TIME 

14400 

SEC 

HMCON 

= 

-4 .6C 

HSUN 

z 

-2.28 

H 

z 

-6.87 

TIME 

= 

1 5000 

SEC 

HMOCN 

-3,93 

HSUN 

= 

-2.67 

H 

= 

-6.61 

TIME 

= 

1 5600 

SEC 

HMOON 

= 

-3.2e 

HSUN 

Z 

-3.05 

H 

z 

-6.31 

TIME 

16200 

SEC 

HMCCN 

-2.59 

ISUN 

z 

-3.39 

H 

Z 

-5.98 

TIME 

1 €800 

SEC 

HMOCN 

-1.92 

HSUN 

Z 

-3.71 

H 

= 

-5.63 

TIME 

s 

1 7400 

SEC 

HMCON 

= 

-1.26 

HSUN 

= 

-4.00 

H 

z 

-5.25 

TIME 

= 

18CC0 

SEC 

HMCCN 

= 

-.61 

HSUN 

z 

-4.25 

H 

z 

-4.86 

TIME 

= 

18600 

SEC 

HMCCN 

= 

.02 

fSUN 

Z 

-4.47 

H 

z 

-4.44 

TIME 

= 

19200 

SEC 

HMOCN 

= 

.63 

HSUN 

z 

-4.65 

H 

z 

-4.02 

TIME 

S 

190  CO 

SEC 

HMOCN 

= 

1.2C 

HSUN 

z 

-4.79 

H 

z 

-3.58 

TIME 

s 

20400 

SEC 

HMCCN 

1.74 

HSUN 

z 

-4.89 

H 

Z 

-3.14 

TIME 

= 

21000 

SEC 

HMCON 

= 

2.24 

fSUN 

z 

-4.94 

h 

Z 

-2.70 

TIME 

= 

21600 

SEC 

HMCON 

2.7C 

HSUN 

z 

-4.96 

H 

s 

-2.26 

TIME 

= 

22200 

SEC 

HMOON 

= 

3.11 

HSUN 

z 

-4.93 

H 

z 

-1.82 
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Table  2.  Surface  Displacement  and  Lunar  and  Solar  Components  of 
Surface  Displacement  for  Station  No.  2 (Continued) 


TIME 

s 

22000 

SEC 

MMOCN 

S 

3.47 

HSCN 

s 

-4.07 

H 

= 

-1.40 

TIME 

= 

23400 

SEC 

HMOON 

= 

3.7  J 

HSUN 

= 

-4.75 

H 

-.90 

TIME 

= 

24000 

SEC 

HMOON 

= 

4.32 

HSUN 

= 

-4.60 

H 

= 

-.50 

TIME 

= 

24E00 

SEC 

HMCCN 

4.21 

HSLN 

= 

-4  .41 

H 

= 

-.20 

TIME 

= 

25200 

SEC 

HMCON 

= 

4.35 

6SUN 

= 

-4.10 

H 

= 

.17 

TIME 

= 

25000 

SEC 

HMCON 

= 

4.42 

HSUN 

-3.91 

H 

= 

.51 

TIME 

= 

264U0 

SEC 

HMCON 

s 

4.42 

HSUN 

= 

-3.61 

H 

= 

.02 

TIME 

= 

270C0 

SEC 

HMCON 

= 

4.37 

HSUN 

-3.27 

H 

= 

1.13 

TIME 

s 

2 ?eoo 

SEC 

HMCCN 

= 

4. 26 

6SUN 

= 

-2.9C 

H 

= 

1.36 

TIME 

= 

20200 

SEC 

HMOON 

= 

4.06 

HSUN 

= 

-2.53 

H 

= 

1.50 

TIME 

s 

20000 

SEC 

HMOON 

= 

3 .05 

HSUN 

= 

-2.00 

H 

= 

1.  77 

TIME 

s 

29400 

SEC 

HMCCN 

= 

3.56 

ISUN 

= 

-1.63 

H 

= 

1.93 

TIME 

= 

30C00 

SEC 

HMCCN 

3.21 

HSUN 

= 

-1 .16 

H 

2.05 

TIME 

= 

30600 

SEC 

HMCCN 

= 

2.01 

HSUN 

= 

-.60 

H 

s 

2.13 

time 

= 

31200 

SEC 

HMCON 

= 

2.37 

HSUN 

r 

-.19 

H 

= 

2.10 

TIME 

= 

310CO 

SEC 

HMCCN 

= 

1.00 

HSUN 

= 

.31 

H 

= 

2.19 

TIME 

s 

32400 

SEC 

HMCON 

= 

1.31 

HSUN 

= 

.02 

H 

s 

2.16 

TIME 

= 

330QO 

SEC 

HMCON 

= 

.77 

HSUN 

r 

1.33 

H 

= 

2.13 

TIME 

s 

33600 

SEC 

HMCCN 

s. 

.17 

(SUN 

r 

1.04 

H 

r 

2.01 

TIME 

=: 

34230 

SEC 

HMCON 

= 

-.46 

HSUN 

= 

2.34 

H 

s 

1.00 

TIME 

= 

340CO 

SEC 

HMOON 

= 

-1.11 

HSUN 

2.03 

H 

= 

1.72 

TIME 

s 

35400 

SEC 

HMCCN 

= 

-1.77 

HSUN 

3.31 

H 

= 

1.53 

TIME 

= 

36000 

SEC 

HMCON 

= 

-2.45 

HSUN 

= 

3.77 

H 

1.32 

TIME 

= 

36600 

SEC 

HMOON 

= 

-3.13 

HSUN 

4.21 

H 

= 

1.00 

TIME 

= 

37200 

SEC 

HMCON 

= 

-3.01 

HSUN 

= 

4.63 

H 

= 

.02 

TIME 

= 

37000 

SEC 

HMOON 

= 

-4.4  1 

HSUN 

= 

5.02 

H 

= 

.53 

TIME 

r 

30403 

SEC 

HMCCN 

= 

-5.14 

6SUN 

= 

5.30 

H 

s 

.23 

TIME 

= 

39000 

SEC 

HMOCN 

s 

-5.7* 

HSUN 

■= 

5.71 

H 

r 

-.30 

TIME 

s 

39600 

SEC 

HMCCN 

= 

-6.4C 

HSUN 

= 

6.33 

H 

= 

-.40 

TIME 

= 

40200 

SEC 

HMCON 

= 

-6.99 

HSUN 

= 

6.26 

H 

-.73 

TIME 

= 

40000 

SEC 

HMCON 

= 

-7.54 

6SUN 

= 

6.40 

H 

= 

-1.06 

TIME 

= 

41400 

SEC 

HMCON 

= 

-0.C5 

HSUN 

s 

6.65 

H 

5 

-1.39 

TIME 

s 

42C0C 

SEC 

HMCON 

= 

-0.51 

HSUN 

= 

6.79 

H 

s 

-1.72 

TIME 

= 

42600 

SEC 

HMCCN 

= 

-0.92 

HSUN 

= 

6.09 

H 

= 

-2.  04 

TIME 

43200 

SEC 

HMCON 

= 

-9.20 

tSUN 

s 

6.94 

H 

= 

-2.35 

TIME 

= 

43003 

SEC 

HMOON 

-9.56 

HSUN 

= 

6.94 

H 

-2.64 

TIME 

s 

44400 

SEC 

HMCON 

= 

-9.02 

HSUN 

= 

6.91 

H 

= 

-2.91 

TIME 

= 

45000 

SEC 

HMCON 

-9.99 

HSUN 

= 

6.03 

H 

r 

-3.16 

TIME 

s 

45600 

SEC 

HMCCN 

= 

-10.09 

>SUN 

r 

6.76 

H 

s 

-3.39 

TIME 

s 

46200 

SEC 

HMOON 

s 

-10.12 

HSUN 

= 

6.54 

H 

= 

-3.50 

TIME 

= 

46000 

SEC 

HMOON 

= 

-10.06 

HSUN 

= 

6.33 

H 

s 

-3.74 

TIME 

47400 

SEC 

HMCCN 

= 

-9.96 

HSUN 

= 

6.09 

H 

= 

-3.07 

TIME 

= 

40000 

SEC 

HMOON 

= 

-9.77 

6SUN 

= 

5.00 

H 

= 

-3.96 

TIME 

= 

40600 

SEC 

HMCON 

= 

-9.5C 

HSUN 

s 

5.49 

H 

s 

-4.02 

TIME 

= 

49230 

SEC 

HMOON 

= 

-9.16 

HSUN 

s 

5.13 

H 

= 

-4.0  3 

TIME 

= 

49000 

SEC 

HMOON 

s 

-0.75 

HSUN 

= 

4.75 

H 

= 

-3.99 

TIME 

r 

50400 

SEC 

HMCCN 

s 

-0.26 

6SUN 

s 

4.34 

H 

= 

-3.91 

TIME 

s 

51C00 

SEC 

HMOON 

= 

-7.7t 

HSUN 

s 

3.91 

H 

-3.79 

time 

s 

51600 

SEC 

HMCCN 

s 

-7.07 

HSUN 

s 

3.46 

H 

= 

-3.62 

TIME 

= 

52200 

SEC 

HMCON 

= 

-6.36 

HSUN 

s 

2.90 

H 

-3.40 

TIME 

= 

52000 

SEC 

HMCON 

= 

-5  .63 

6 SUN 

s 

2.56 

H 

= 

-3.13 

TIME 

s 

5 3400 

SEC 

HMCON 

r 

-4.01 

HSUN 

= 

2.00 

H 

= 

-2.01 

TIME 

s 

54CC0 

SEC 

HMOON 

s 

-3.95 

HSUN 

= 

1.49 

H 

s 

-2.45 

TIME 

= 

54600 

SEC 

HMCCN 

= 

-3.03 

HSUN 

= 

.99 

H 

s 

-2.05 

TIME 

= 

55200 

SEC 

HMCON 

= 

-2.07 

HSUN 

s 

.40 

H 

= 

-1.60 

J 


Table  2.  Surface  Displacement  and  Lunar  and  Solar  Components  of 
Surface  Displacement  for  Station  No.  2 (Continued) 


TIME 

s 

55800 

SEC 

hmcon 

= 

-1.07 

HSUN 

-.03 

H 

X 

-1.10 

TIHE 

= 

56400 

SEC 

HMCCN 

= 

-.04 

hSUN 

-.52 

H 

X 

-.57 

TIME 

s 

570GO 

SEC 

HMCON 

= 

1.01 

HSON 

-1.01 

H 

X 

.01 

TIME 

X 

57600 

SEC 

HMCON 

= 

2.05 

HSUN 

= 

-1.48 

H 

X 

.61 

TIME 

S 

56200 

SEC 

HMCON 

3.1  * 

HSUN 

= 

-1.93 

H 

X 

1.25 

TIME 

s 

5 8803 

SEC 

HMCCN 

= 

4.2* 

HSLN 

= 

-2.36 

H 

X 

1.92 

TIME 

= 

59400 

SEC 

HMCON 

= 

5.39 

fSUN 

X 

-2.77 

H 

X 

2. 62 

TIME 

= 

600C0 

SEC 

HMCON 

6.4  1 

HSUN 

X 

-3.15 

H 

X 

3.34 

TIME 

s 

606  00 

SEC 

hmcon 

= 

7.57 

HSUN 

X 

-3.50 

H 

X 

4.07 

time 

= 

6 1200 

SEC 

HMCON 

X 

8.63 

HSUN 

X 

-3.81 

H 

X 

4.82 

TIME 

s 

618Q0 

SEC 

HMCCN 

9 .67 

ESUN 

X 

-4.09 

H 

X 

5.58 

TIME 

s 

6240Q 

SEC 

HMCON 

= 

10.6  ( 

HSUN 

X 

-4.34 

H 

X 

6.34 

TIME 

= 

63000 

SEC 

HMCON 

X 

11.64 

HSUN 

X 

-4.54 

H 

X 

7.10 

TIME 

= 

63600 

SEC 

HMCCN 

X 

12.57 

HSUN 

X 

-4.71 

H 

X 

7.86 

TIME 

64200 

SEC 

HMCCN 

X 

13  .44 

ESUN 

X 

-4.83 

H 

X 

8.61 

TIME 

= 

64800 

SEC 

HMCCN 

X 

14. 2( 

HSUN 

X 

-4  .91 

H 

X 

9.34 

TIME 

= 

65400 

SEC 

HMCON 

X 

15.01 

HSUN 

X 

-4.95 

H 

X 

10.06 

TIME 

s 

66000 

SEC 

HMCON 

X 

15.70 

HSUN 

X 

-4.95 

H 

X 

10.75 

TIME 

s 

66600 

SEC 

HMCON 

X 

16.32 

ESUN 

X 

-4.91 

H 

X 

11.41 

time 

= 

67200 

SEC 

HMCON 

X 

16.86 

HSUN 

X 

-4.82 

E 

X 

12.04 

TIME 

X 

67830 

SEC 

HMCON 

X 

17.33 

HSLN 

X 

-4.70 

H 

X 

12.63 

TIME 

X 

68400 

SEC 

HMCCN 

X 

17.71 

ESUN 

X 

-4.53 

H 

X 

13.18 

TIME 

X 

65000 

SEC 

HMCON 

X 

18.0  1 

HSUN 

X 

-4.33 

H 

X 

13.68 

time 

X 

69600 

SEC 

HMCON 

X 

18.23 

HSUN 

X 

-4.09 

H 

X 

14.14 

TIME 

X 

702C0 

SEC 

HMCCN 

X 

18.35 

HSUN 

X 

-3.81 

H 

X 

14.54 

TIME 

X 

7 0000 

SEC 

HMCON 

X 

18.39 

HSUN 

X 

-3.51 

H 

X 

14.89 

TIME 

X 

71400 

SEC 

HMCCN 

X 

18.35 

HSUN 

X 

-3.17 

H 

X 

15.17 

TIME 

X 

72000 

SEC 

HMCON 

X 

18.21 

HSUN 

X 

-2.81 

H 

X 

15.40 

TIME 

X 

7 2600 

SEC 

HMCON 

X 

17.99 

HSUN 

X 

-2.42 

H 

X 

15.56 

TIME 

X 

73200 

SEC 

HMCON 

X 

17.68 

ESUN 

X 

-2.02 

H 

X 

15.66 

TIME 

s 

73800 

SEC 

HMOCN 

X 

17.2  1 

HSUN 

X 

-1.55 

H 

X 

15.69 

TIME 

X 

74400 

SEC 

HMCON 

X 

16.81 

HSUN 

X 

-1.15 

H 

X 

15.65 

TIME 

X 

75CC0 

SEC 

HMCON 

X 

16.25 

HSUN 

X 

-.71 

H 

X 

15.55 

TIME 

= 

75600 

SEC 

HMOCN 

X 

15.62 

ESUN 

X 

-.25 

H 

X 

15.37 

time 

s 

76200 

SEC 

HMCON 

X 

14. 9< 

HSUN 

X 

.21 

H 

X 

15.13 

TIME 

X 

76800 

SEC 

HMCON 

X 

14. If 

HSUN 

X 

.67 

E 

X 

14.83 

time 

X 

77400 

SEC 

HMCCN 

X 

13.33 

HSUN 

X 

1.12 

H 

X 

14.45 

TIME 

= 

78C00 

SEC 

HMCON 

X 

12.44 

ESUN 

X 

1.57 

H 

X 

14.01 

TIME 

= 

786CQ 

SEC 

HMOCN 

X 

11. 5C 

HSUN 

X 

2.01 

E 

X 

13.51 

TIME 

s 

75200 

SEC 

HMOCN 

X 

10.52 

HSUN 

X 

2.43 

h 

X 

12.95 

time 

X 

758C0 

SEC 

HMCON 

X 

9.50 

HSUN 

X 

2.83 

H 

X 

12.34 

TIME 

X 

80400 

SEC 

HMCCN 

X 

8.45 

ESUN 

X 

3.22 

H 

X 

11.67 

TIME 

X 

81000 

SEC 

HMCON 

X 

7.37 

HSUN 

X 

3.58 

H 

X 

10.95 

TIME 

X 

81600 

SEC 

HMOCN 

X 

6.27 

HSUN 

X 

3.91 

H 

X 

10.18 

TIME 

X 

82200 

SEC 

HMCON 

X 

5. If 

HSUN 

X 

4.21 

H 

X 

9.37 

TIME 

X 

8 2800 

SEC 

HMOON 

X 

4.04 

ESUN 

X 

4.48 

H 

X 

8.52 

TIME 

X 

83400 

SEC 

HMOCN 

X 

2.92 

HSUN 

X 

4.72 

H 

X 

7.64 

TIME 

X 

84000 

SEC 

HMOON 

X 

1.82 

HSUN 

X 

4.92 

H 

X 

6.74 

time 

X 

84600 

SEC 

HMOON 

X 

.72 

HSUN 

X 

5.09 

H 

X 

5.81 

TIME 

X 

85200 

SEC 

HMCCN 

X 

-.35 

ESUN 

X 

5.21 

H 

X 

4.86 

TIME 

X 

85800 

SEC 

HMOCN 

X 

-1.35 

HSUN 

X 

5.30 

H 

X 

3.90 

25 


Table  3.  Surface  Displacement  and  Lunar  and  Solar  Components  of 
Surface  Displacement  for  Station  No.  3 


(Station  longitude  = 0 deg.  station  latitude  = 60  deg.  HMOON,  HSUN,  and  H are 
in  centimeters.  Time  is  listed  in  terms  of  seconds  from  Midnight  UT  of  Day  88  of 
Year  1977.) 


TIME 

s 

a 

SEC 

HMOON 

-5.75 

HSUN 

s 

-1.95 

H 

s 

-7.69 

TIME 

= 

600 

SEC 

HMCCN 

= 

-6.11 

HSLN 

= 

-1 .95 

H 

r 

-8.12 

TIME 

= 

1200 

SEC 

HMCCN 

= 

-6.59 

ISON 

= 

-1.96 

H 

= 

-8.55 

TIME 

1600 

SEC 

HMOON 

-6.96 

HSUN 

-1.98 

H 

-8.96 

TIME 

2i«Q0 

SEC 

HMOON 

= 

-7.35 

HSUN 

= 

-2.02 

H 

= 

-9.37 

TIME 

= 

3000 

SEC 

HMCCN 

= 

-7  .E« 

HSLN 

— 

-2.07 

H 

s 

-9.76 

TIME 

- 

3600 

SEC 

HMOON 

-8.02 

► SUN 

-2.13 

H 

X 

-10,14 

TIME 

= 

^200 

SEC 

HMCCN 

= 

-8.31 

HSUN 

-2.20 

H 

r 

-10.51 

TIME 

= 

4600 

SEC 

HMOON 

= 

-8.5  « 

hSLN 

-2.28 

H 

= 

-1C. 87 

TIME 

5400 

SEC 

HMCCN 

= 

-8.84 

HSUN 

-2.37 

H 

-11.21 

TIME 

= 

6000 

SEC 

HMOON 

■= 

-9.06 

HSUN 

= 

-2.47 

H 

-11.53 

TIME 

s 

6600 

SEC 

HMCON 

-9.25 

HSUN 

= 

-2  .56 

H 

= 

-11.84 

TIME 

= 

7200 

SEC 

HMCCN 

= 

-9.43 

HSLN 

-2  .70 

H 

s 

-12.13 

TIME 

s 

7600 

SEC 

HMCON 

-9.57 

>SUN 

-2.82 

H 

= 

-12.40 

TIME 

= 

6400 

SEC 

HMOON 

-9.6  ♦ 

HSUN 

= 

-2.95 

H 

= 

-12. 65 

TIME 

= 

9000 

SEC 

HMCON 

= 

-9.7  < 

HSUN 

-3.09 

H 

-12.88 

TIME 

9600 

SEC 

HMCCN 

-9.87 

HSUN 

= 

-3.22 

H 

= 

-13.09 

TIME 

= 

1 0200 

SEC 

HMCCN 

- 

-9.92 

ISUN 

-3.36 

H 

= 

-13.28 

TIME 

1 06  30 

SEC 

HMOON 

-9.95 

HSUN 

= 

-3.50 

H 

-13.45 

TIME 

= 

11430 

SEC 

HMCON 

-9  .96 

HSUN 

s 

-3.64 

H 

= 

-13.61 

TIME 

= 

12000 

SEC 

HMCON 

= 

-9.96 

HSLN 

= 

-3.78 

H 

-13.74 

TIME 

12600 

SEC 

HMCON 

-9  .93 

HSUN 

-3.91 

H 

= 

-13.65 

TIME 

= 

13200 

SEC 

HMCCN 

= 

-9.9( 

HSUN 

- 

-4.05 

H 

-13.94 

TIME 

' = 

136  30 

SEC 

HMOON 

= 

-9.84 

HSUN 

s 

-4.17 

H 

= 

-14. 02 

TIME 

14400 

SEC 

HMCON 

-9.76 

HSUN 

= 

-4.29 

H 

-14. C7 

TIME 

= 

15000 

SEC 

HMCCN 

= 

-9.7C 

HSUN 

= 

-4.41 

H 

-= 

-14.11 

TIME 

= 

156C0 

SEC 

HMOON 

-9,61 

HSUN 

-4.51 

H 

= 

-14.13 

TIME 

= 

162C0 

SEC 

HMCON 

= 

-9.52 

HSUN 

= 

-4.61 

H 

= 

-14.13 

TIME 

16603 

SEC 

HMCON 

= 

-9  .42 

fSUN 

s 

-4.69 

H 

s 

-14.12 

TIME 

= 

1 7400 

SEC 

HMOON 

= 

-9.32 

HSUN 

= 

-4.77 

H 

= 

-14. C9 

TIME 

s 

1 6000 

SEC 

HMOON 

:: 

-9.22 

HSUN 

s 

-4.83 

H 

s 

-14.05 

TIME 

s 

166C0 

SEC 

HMCCN 

= 

-9.11 

HSUN 

= 

-4.88 

H 

= 

-14.03 

TIME 

= 

19200 

SEC 

HMCON 

= 

-9.01 

HSUN 

s 

-4.92 

H 

-13.93 

TIME 

= 

19800 

SEC 

HMOON 

= 

-8.91 

HSUN 

s 

-4.95 

H 

= 

-13.85 

TIME 

= 

20400 

SEC 

HMCON 

= 

-8.81 

HSUN 

= 

-4.96 

H 

= 

-13.77 

TIME 

= 

210CQ 

SEC 

HMOON 

-8,72 

HSUN 

2: 

-4.96 

H 

-13.67 

TIME 

s 

21600 

SEC 

HMCON 

-8.63 

ISUN 

= 

-4  .94 

H 

s 

-13.57 

TIME 

= 

22200 

SEC 

HMOON 

= 

-8.55 

HSUN 

-4.91 

H 

£ 

-13.46 

Table  3.  Surface  Displacement  and  Lunar  and  Solar  Components  of 
Surface  Displacement  for  Station  No.  3 (Continued) 


TIHE 

s 

2ZSQ0 

SEC 

HMCCN 

r 

-8.4  1 

hSUN 

= 

-4.87 

H 

s 

-13.35 

TIME 

= 

Z34a0 

SEC 

HMOON 

= 

-8.42 

HStIN 

= 

-4.81 

H 

3 

-12.23 

TIME 

24CQ0 

EEC 

HNCON 

s 

-8  .37 

fSUN 

= 

-4.73 

H 

= 

-13.11 

TIME 

= 

24eoc 

SEC 

HMCCN 

= 

-8. 24 

HSUN 

= 

-4.65 

H 

= 

-12.99 

TIME 

= 

25200 

SEC 

HMCCN 

= 

-8  .31 

HSUN 

= 

-4 .55 

H 

3 

-12.86 

TIME 

25000 

EEC 

HMCCN 

r 

-8.30 

HSUN 

= 

-4.44 

H 

= 

-12.74 

TIHE 

= 

26400 

EEC 

HMCCN 

= 

-8.3C 

HSUN 

= 

-4.32 

H 

= 

-12.61 

TIME 

= 

27000 

SEC 

HMCCN 

s 

-8.31 

hSUN 

= 

-4.18 

H 

S 

-12.49 

TIME 

= 

27600 

SEC 

HMCCN 

= 

-8.33 

HSUN 

= 

-4.04 

H 

= 

-12.37 

TIME 

= 

20200 

SEC 

HMCCN 

= 

-8.37 

hSUN 

= 

-3.89 

H 

= 

-12.25 

TIME 

28800 

SEC 

HMCCN 

= 

-8.41 

HSUN 

s 

-3.72 

H 

3 

-12.14 

TIKE 

= 

29400 

SEC 

HMCCN 

= 

-8.47 

hSUN 

r 

-3.56 

H 

S 

-12.03 

TIME 

= 

30CC0 

SEC 

HMCCN 

= 

-8.54 

HSUN 

= 

-3.38 

H 

= 

-11.92 

TIME 

= 

30600 

SEC 

HMCCN 

s 

-8.62 

HSUN 

= 

-J.2C 

H 

= 

-11.82 

TIME 

s 

31200 

SEC 

HMCCN 

= 

-8.70 

hSUN 

r 

-3.02 

h 

3 

-11.72 

TIME 

= 

31800 

SEC 

HMCCN 

= 

-8.80 

HSUN 

= 

-2.84 

H 

3 

-11.63 

TIME 

= 

32400 

SEC 

HMCCN 

= 

-e  .9( 

hSUN 

= 

-2.65 

H 

3 

-11.55 

TIME 

s 

33000 

SEC 

HMCCN 

= 

-9.0C 

HSUN 

s 

-2.46 

H 

S 

-11.46 

TIME 

s 

33600 

SEC 

HMCCN 

s 

-9.11 

HSUN 

s 

-2.28 

h 

3 

-11.39 

TIME 

= 

34200 

SEC 

HMCCN 

s 

-9.22 

HSUN 

= 

-2.10 

H 

3 

-11.32 

TIME 

= 

34800 

SEC 

HMCCN 

= 

-9.32 

hSUN 

s 

-1.92 

H 

3 

-11.25 

TIME 

= 

35400 

SEC 

HMCCN 

■= 

-9.44 

HSUN 

= 

-1.75 

H 

3 

-11.19 

TIME 

s 

36000 

SEC 

HMCCN 

= 

-9.55 

HSUN 

= 

-1.58 

h 

3 

-11.13 

TIME 

s 

36600 

SEC 

HMCCN 

= 

-9.65 

HSUN 

3 

-1.42 

H 

3 

-11.07 

TIME 

s 

37200 

SEC 

HMCCN 

-9.74 

hSUN 

= 

-1.27 

H 

3 

-11. C2 

TIME 

= 

37800 

SEC 

HMCCN 

= 

-9.8  2 

HSUN 

3 

-1.13 

H 

3 

-10.96 

TIME 

r 

38400 

SEC 

HMOON 

s 

-9.91 

HSUN 

s 

-l.OC 

h 

3 

-10.91 

TIME 

= 

39000 

SEC 

H.'ICCN 

s 

-9.98 

HSUN 

3 

-.89 

H 

3 

-10.86 

TIME 

= 

39600 

SEC 

HMCCN 

= 

-10.03 

hSUN 

= 

-.76 

H 

3 

-10.81 

TIME 

= 

4C2C0 

SEC 

HMCCN 

= 

-10.07 

HSUN 

s 

-.69 

H 

3 

-1C. 76 

TIME 

= 

40800 

SEC 

HMCCN 

-10 .0  • 

HSUN 

= 

-.61 

h 

3 

-10.70 

TIME 

5 

41400 

SEC 

HMCCN 

= 

-10. 1C 

HSUN 

= 

-.55 

H 

3 

-10.64 

TIME 

= 

42000 

SEC 

HMCCN 

: 

-10.08 

hSUN 

= 

-.50 

H 

= 

-10.58 

time 

s 

42600 

SEC 

HMOON 

-10.09 

HSUN 

= 

-.46 

H 

3 

-10. 51 

TIME 

= 

43200 

SEC 

HMCCN 

r 

-9.99 

HSUN 

s 

-.44 

H 

3 

-10.44 

TIME 

= 

43800 

SEC 

HMCCN 

= 

-9.92 

HSUN 

= 

-.44 

H 

3 

-10.36 

TIME 

= 

44400 

SEC 

HMCCN 

= 

-9.81 

HSUN 

= 

- .45 

H 

3 

-10.27 

TIME 

= 

45000 

SEC 

HMCCN 

= 

-9.69 

HSUN 

s 

-.48 

H 

3 

-10.17 

TIME 

s 

456C0 

SEC 

HMCCN 

= 

-9.54 

HSUN 

-.52 

H 

3 

-10.06 

TIME 

= 

46200 

SEC 

HMCCN 

= 

-9 .3e 

hSUN 

3 

-.58 

H 

3 

-9.95 

TIME 

s 

46800 

SEC 

HMCCN 

= 

-9. If 

HSUN 

3 

-.66 

H 

3 

-9.82 

TIME 

= 

47400 

SEC 

HMCCN 

-8.94 

HSUN 

= 

-.74 

h 

3 

-9.68 

TIME 

48000 

SEC 

HMCCN 

= 

-8.69 

HSUN 

3 

-.84 

H 

3 

-9.53 

TIME 

= 

48600 

SEC 

HMCCN 

= 

-8.41 

ISUN 

3 

-.9b 

H 

3 

-9.37 

TIME 

= 

49200 

SEC 

HMCCN 

-8.1  3 

HSUN 

3 

-1.08 

H 

3 

-9.19 

TIME 

= 

49800 

SEC 

HMCCN 

= 

-7.7< 

HSUN 

s 

-1.22 

H 

3 

-9. Cl 

TIME 

= 

50400 

SEC 

HMCCN 

= 

-7.44 

HSUN 

3 

-1.36 

H 

3 

-8.81 

TIME 

= 

51CC0 

SEC 

HMCCN 

s 

-7.07 

hSUN 

3 

-1.52 

H 

3 

-8.59 

TIME 

= 

51603 

SEC 

HMCCN 

= 

-6.6t 

HSUN 

3 

-1.6  f 

H 

3 

-8.36 

TIME 

s 

52200 

SEC 

HMCCN 

-6.27 

HSUN 

3 

-1.85 

H 

3 

-8.12 

TIME 

= 

52003 

SEC 

HMCCN 

= 

-5.84 

HSUN 

= 

-2.03 

H 

3 

-7.87 

TIME 

= 

53400 

SEC 

HMOON 

r 

-5.40 

hSUN 

3 

-2.21 

H 

3 

-7.61 

TIME 

= 

54000 

SEC 

HMOON 

-4.94 

HSUN 

3 

-2.39 

H 

3 

-7.33 

time 

s 

54600 

SEC 

HMOON 

= 

-4  .46 

HSUN 

= 

-2.58 

H 

3 

-7.04 

time 

= 

55200 

SEC 

HMCCN 

= 

-3.98 

HSUN 

= 

-2.76 

H 

3 

-6. 7'* 
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Table  3.  Surface  Displacement  and  Lunar  and  Solar  Components  of 
Surface  Displacement  for  Station  No.  3 (Continued) 


TIME 

55«Ca 

SEC 

HMCCN 

-3.49 

ESUN 

= 

-2.95 

H 

= 

-6.43 

TIME 

56^400 

SEC 

HMC0.9 

= 

-2.95 

HSJN 

s 

-3.13 

H 

= 

-6.12 

TIME 

= 

57GOO 

sec 

HMCCN 

= 

-2.4  f 

HSUN 

= 

-3.31 

H 

s 

-5.79 

TIME 

= 

5 TECO 

sec 

HMCCM 

= 

-1.97 

HSCN 

-3.45 

I- 

S 

-5.46 

TIME 

= 

5 ezoa 

sec 

HMCO.'J 

= 

-1.47 

HSUN 

= 

- 3 . 6(3 

H 

s 

-5.13 

TIME 

s 

5 flsao 

sec 

HriOON 

= 

-.96 

HSUN 

= 

-3.82 

6 

r 

-4.79 

TIME 

= 

59<4Ca 

sec 

HMCCM 

= 

- .47 

HSUN 

= 

-3.98 

H 

s 

-4. 4m 

TIME 

= 

eocoa 

SEC 

HMCCN 

= 

.03 

iSUN 

= 

-4.12 

H 

= 

-4.10 

TIME 

= 

6 06  00 

sec 

HMOON 

.5 : 

HSUN 

= 

-4.26 

H 

-3.76 

TIME 

= 

61tC0 

SEC 

HMOCN 

= 

.97 

HSUN 

= 

-4.39 

H 

= 

-3.41 

TIME 

= 

6 18C3 

SEC 

HMCCN 

= 

1.43 

HSUN 

s 

-4.53 

H 

-3.08 

TIME 

= 

62<4CQ 

SEC 

HMCCN 

= 

1 .86 

6SUN 

= 

-4.61 

H 

= 

-2.  74 

TIME 

63GC0 

sec 

HMCCN 

= 

2.2( 

HSUN 

-4  .73 

H 

= 

-2.42 

TIME 

= 

636  CO 

SEC 

HMOON 

= 

2.67 

HSUN 

= 

-4.76 

H 

= 

-2. 10 

TIME 

= 

64203 

sec 

HMCCN 

= 

3.04 

HSUN 

= 

-M  . 84 

H 

= 

-1.80 

TIME 

= 

64600 

sec 

H.MCON 

= 

3.35 

ESUN 

= 

-4 . 89 

H 

= 

-1.50 

time 

= 

65400 

SEC 

HMCCN 

= 

3.71 

HSUN 

-4.93 

H 

r 

-1.22 

time 

= 

66003 

SEC 

HMOON 

3.95 

HSUN 

= 

-4.95 

H 

X 

-.96 

TIKE 

r 

66600 

SEC 

HMCCN 

- 

4.25 

HSUN 

= 

-4.gt 

H 

X 

-.71 

TIME 

= 

672C0 

sec 

HMCCN 

r 

4.47 

6 SUN 

-4.95 

H 

= 

-.48 

TIME 

= 

67603 

SEC 

HMCCN 

= 

4.66 

H jUN 

= 

-4.93 

H 

= 

-.27 

TIME 

s 

66403 

sec 

HMCCN 

= 

H . 8 2 

HSUN 

r 

-4  .9C 

h 

X 

-.06 

TIME 

s 

69CC0 

SEC 

H ICON 

= 

4.94 

HSUN 

= 

-4.86 

H 

X 

.08 

TIME 

= 

69630 

SEC 

HMCON 

5 .C2 

ISUN 

-4.80 

H 

X 

.22 

TIME 

7Q2C0 

SEC 

HMCCN 

= 

5.36 

HSUN 

= 

-4.73 

H 

X 

.34 

TIME 

= 

706C0 

SEC 

HMCCN 

= 

5. 0 ; 

HSUN 

= 

-4.65 

H 

= 

.43 

TIME 

71400 

SEC 

HMCCN 

5.04 

HSUN 

-4.56 

h 

= 

.49 

TIME 

= 

720  00 

SEC 

HMCCN 

4.98 

iSUN 

= 

-4.46 

H 

X 

.52 

TIME 

•= 

7260U 

SEC 

HMCON 

= 

4.83 

HSUN 

= 

-4.35 

H 

= 

.53 

TIME 

73203 

SEC 

HMCCN 

= 

4 .74 

HSUN 

s 

-4.23 

H 

X 

.50 

TIME 

= 

73603 

sec 

HMOON 

= 

4. 56 

HSUN 

= 

-4.11 

H 

X 

.45 

TIME 

s 

744  00 

SEC 

HMCCN 

= 

4.3! 

ISUN 

= 

-3.96 

H 

X 

.37 

TIME 

= 

75003 

SEC 

HMCC  i 

= 

4.i: 

HSUN 

z. 

-3.65 

H 

X 

.26 

TIME 

756C0 

SEC 

HMCCN 

= 

3.8i 

HSUN 

= 

-3.72 

h 

X 

.12 

time 

= 

76233 

sec 

HMCCN 

= 

3.5! 

HSUN 

-3.58 

H 

X 

-.  Cb 

TIME 

= 

76600 

SEC 

HMOON 

3.1? 

ISUN 

= 

-3  .h5 

H 

X 

-.25 

time 

r 

77400 

SEC 

HMCCN 

r 

2.82 

HSUN 

z 

-3.31 

H 

X 

-.48 

time 

= 

7eC30 

sec 

HMCON 

= 

2.44 

HSUN 

= 

-3.17 

6 

X 

-.74 

TIME 

7 6633 

SEC 

HMCCN 

= 

2.02 

HSUN 

= 

-3.C4 

H 

X 

-1.02 

TIME 

= 

7 92  CO 

SEC 

HMCCN 

= 

1.  is 

E SUN 

-2.91 

H 

= 

-1.33 

TIME 

= 

7 58C0 

SEC 

HMCCN 

= 

1.14 

HSUN 

= 

-2.79 

H 

= 

-1 . €o 

TIME 

= 

8 3403 

SEC 

HMCCN 

s 

. >6 

HSUN 

= 

-2.6  6 

H 

X 

-2.01 

TIME 

r 

81030 

SEC 

HMCCN 

s 

.16 

ISUN 

= 

-2.57 

h 

-2.39 

TIME 

s 

816uJ 

SEC 

HMtoN 

= 

-.32 

HSUN 

r 

-2.47 

H 

r 

-2.  78 

TIME 

s 

8 2200 

SEC 

HMCCN 

= 

-.63 

HSUN 

= 

-2.37 

H 

= 

-3.20 

TIME 

= 

82633 

SfcC 

HMOCN 

= 

-1.34 

HSUN 

= 

-2.29 

H 

X 

-3.63 

TIME 

83400 

SEC 

HMCCN 

r 

-1.66 

6 SUN 

= 

-2.22 

H 

= 

-4.08 

time 

= 

84C3j 

SEC 

HM  OON 

= 

-2.36 

HSUN 

= 

-2.16 

H 

= 

-4. 5m 

time 

= 

846  3C 

SEC 

HMCCN 

= 

-2.9: 

HSUN 

s 

-2.11 

H 

X 

-5.  :i 

time 

= 

85200 

SEC 

HMCCN 

= 

-3.42 

HSUN 

= 

-?.07 

H 

z 

-5.49 

time 

3 

856C3 

SEC 

HMCCN 

r 

-3.93 

ISUN 

r 

-2.04 

H 

X 

-5.57 

28 


Table  4.  Sample  Data  from  Two-Hand  Computations  for  Computer 

Program  Checkout 


STATION  # 

1 

2 

STN  LONG  (deg) 

0 

0 

STN  LAT  (deg) 

0 

30 

UNIVERSAL  TIME  (sec) 

OF  DAY  88  IN  YEAR  1977 

57600 

1200 

XSUN  (km)  = Xg 

77220921 

-148806042 

YSUN  = y's 

-127563246 

9783868 

ZSUN  = Zg 

9143321 

8480080 

XMOON  (km)  = x'^ 

229338 

65180 

YMOON  = y'l 

300370 

-372772 

ZMOON  = z'l 

103334 

113461 

RSUN  (km) 

149395625 

149368250 

^SUN  (deg) 

3.5088129 

3.254602 

XSUN  (deg) 

-58.3934085 

1 76.656074 

RMOON  (km) 

391786 

395071 

i^MOON  (deg) 

15.292773 

16.689852 

XMOON  (deg) 

53.055410 

-79.6642185 

cos  y SUN 

0.523101438 

-0.834769991 

7 SUN  (deg) 

58.459 

146.591896 

(cos  7 SUN) 

-0.08954733 

0.545261408 

cos  7 MOON 

0.579759898 

0.292429179 

7 MOON  (deg) 

54.566343 

72.9965565 

Pj  (cos  7 MOON) 

0.004182309 

-0.371727763 

HSUN  (cm) 

-0.9 

5.4 

HMOON  (cm) 

0.085 

-7.4 

{ 
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The  following  are  explanations  for  some  of  the  quantities  listed  in  Table  4. 


XSUN,  YSUN,  ZSUN  and  XMOON,  YMOON,  ZMOON  are  the  earth-fixed  Cartesian 
coordinates  of  the  “fictitious”  sun  and  moon  at  the  specified  Universal  time.  These 
were  produced  by  the  computer  algorithm  while  excecuting  those  computations 
which  led  to  the  results  listed  in  Tables  1 and  2.  They  were  made  the  starting 
point  of  the  two  hand  cases  without  being  recomputed  by  hand,  from  the  sun  and 
moon  ephemerides.  Various  checks  were  however  made  to  assure  correctness,  such  as 
comparison  with  the  positions  listed  in  the  American  Ephemeris  and  Nautical 
Almanac  as  well  as  tests  concerning  the  direction  and  rate  of  rotation  of  the 
fictitious  sun  and  moon  with  respect  to  inertial  space  and  with  respect  to  the  earth. 

RSUN,  SUN,  and  X SUN  are  the  geocentric  distance,  latitude  and  longitude  of  the 
fictitious  sun.  RMOON,  ipMOON,  XMOON  are  the  corresponding  quantities  for  the 
fictitious  moon. 

7 SUN  is  the  spherical  angle  between  the  station  and  the  fictitious  sun.  7 MOON  is 
the  corresponding  quantity  for  the  fictitious  moon. 

HSUN  and  HMOON  are  the  surface  displacements  caused  by  sun  and  moon,  at  the 
specified  value  of  Universal  time. 


FUTURE  REQUIREMENTS 

The  computer  program  specified  in  this  report  for  the  tidal  surface  elevation  is 
regarded  adequate  for  use  with  the  present  version  of  CELEST.  But  now  already  a 
number  of  desirable  improvements  may  be  anticipated. 

The  text  following  Equation  (15)  relates  that  the  tide  potential  and  thus  the 
surface  displacement  may  be  truncated  after  the  term  containing  hj.  But  as 
Appendix  A will  demonstrate,  even  this  latter  term  may  safely  be  deleted  unless 
millimeter  accuracy  is  desired.  On  the  other  hand,  Reference  6,  Page  7 relates  that 
the  orbit  integration  for  geodetic  satellites  already  makes  use  of  the  dynamical  Love  j 

numbers  up  to  kj  and  also  resorts  to  surface  harmonic  expansions  for  the  Love 
numbers  in  order  to  account  for  the  details  of  the  gravitational  influence  of  the  ' 

tidal  bulge  on  the  satellite  orbit.  i 


The  present  version  of  TERRA  considers  the  dynamical  Love  numbers  program 
input.  But  the  author  expects  that  in  the  not  too  distant  future  these  numbers  may 
be  made  bias  parameters  to  be  “measured”  by  the  least  squares  procedure  involved 
in  the  TERRA  runs.  At  present,  CELEST  also  treats  the  Love  numbers  hj  and 
perhaps  hj  as  input  quantities.  Consequently,  it  should  be  interesting  to  expand  h^ 
into  surface  harmonics  and  to  experiment  with  the  thus  modified  algorithm.  And, 
secondly,  in  the  liglit  of  Appendix  C,  one  may  attempt  to  establish  a valid 
relationship  amongst  the  dynamical  Love  numbers  resulting  from  the  geodetic 
parameter  solutions  and  the  yet  unmeasured  expansion  coefficients  for  the 
geometrical  Love  numbers  hj  and  lij.  The  author  feels  certain  that  the  necessary 
relationship  h^  = ) can  be  approximated  by  a rather  straightforward  extension 

of  the  elementary  theory  involving  the  geometrical  and  dynamical  Love  number 
concepts.  It  is  believed  that  this  extension  will  involve  essentially  a mass  density 
model  which  deviates  from  the  customary  assumption  of  homogeneity.  The  author 
plans  to  experiment  with  a three-layer  model  of  the  tidal  earth.  Also,  a review  is 
being  made  of  the  latest  geophysical  research  relevant  to  the  solid  earth  tide,  with 
particular  attention  to  work  published  by  P.  Melchior  (note  the  summary  of  the 
subject  in  Reference  10,  Volume  3). 


REFERENCES 


1.  J.  W.  O’Toole,  Functional  Description  of  Satellite  Ephemerides  Program 
CELEST,  Naval  Surface  Weapons  Center,  Dahlgren  Laboratory  Technical  Report 
NSWe/DL  TR-2849,  Dahlgren,  VA,  July  1976. 

2.  J.  W.  O’Toole,  CELEST  Computer  Program  for  Computing  Satellite  Orbits, 
Naval  Surface  Weapons  Center,  Dahlgren  Laboratory  Technical  Report  NSWC/DL 
TR-3565,  Dahlgren,  VA,  October  1976. 

3.  W.  Groeger  and  R.  Manrique,  Ocean  Tides  in  the  Equations  of  Motion  of  the 
TERRA  System  of  Computer  Programs  for  Satellite  Geodesy,  Naval  Surface 
Weapons  Center,  Dahlgren  Laboratory  Technical  Report  NSWC/DL  TR-3535j 
Dahlgren,  VA,  December  1976. 

4.  W.  Groeger,  A Gravitational  Potential  for  Atmospheric  Earth  Tides  Caused  by 
the  Sun,  Naval  Surface  Weapons  Center,  Dahlgren  Laboratory  Technical  Report 
NSWC/DL  TR-3485,  Dahlgren,  VA,  October  1976.  ' 


31 


5.  R.  Manrique  and  W.  Groeger,  A Gravitational  Potential  for  Atmospheric  Earth 
Tides  Caused  by  the  Moon,  Naval  Surface  Weapons  Center,  Dahlgren  Laboratory 
Technical  Report  NSWC/DL  TR-3638,  Dahlgren,  VA,  November  1976. 

6.  W.  J.  Groeger  and  R.  B.  Manrique,  Solid  Earth  Tides  in  the  Equations  of 
Motion  of  the  TERRA  System  of  Computer  Programs  for  Satellite  Geodesy, 
Naval  Surface  Weapons  Center,  Dahlgren  Laboratory  Technical  Report  NSWC/DL 
TR-3560,  Dahlgren,  VA,  February  1977. 

7.  H.  Jeffreys,  The  Earth,  Cambridge  University  Press,  Cambridge,  5“’  Edition, 
1970. 

8.  Z.  Kopal,  Physics  and  Astronomy  of  the  Moon,  Academic  Press,  New  York  and 
London,  1962. 

9.  K.  Lambeck  and  A.  Cazenave,  Fluid  Tidal  Effects  on  Satellite  Orbit  and  Other 
Temporal  Variations  in  the  Geopotential,  Bulletin  Groupc  Rcch.  Geod.  Spatiale, 
7,  20,  January  1973. 

10.  P.  Melchior,  Physique  et  Dynamique  Planetaires,  3 Volumes,  Vander-Editeur, 
Louvain  and  Bruxelles,  1971  - 1972. 

11.  W.  H.  Munk  and  J.  F.  MacDonald,  The  Rotation  of  the  Earth,  Cambridge 
University  Press.  Cambridge,  1960. 

12.  P.  Melchior,  Past  and  Future  of  Research  Methods  in  Problems  of  the  Earth's 
Rotation,  Communications,  Series  B,  No.  71,  Observatoire  Royal  de  Belgique, 
Reprinted  from  Rotation  of  the  Earth  lAU  Symposium  48,  Morioka.  Japan, 
1972. 

13.  P.  Melchior,  Earth  Tides  and  Polar  Motions,  Communications,  Series  B,  No.  76, 
Observatoire  Royal  de  Belgique,  Reprinted  from  Tectonophysics,  Vol.  13,  1972. 

14.  P.  Melchior  and  B.  Georis,  Earth  Tides,  Precession-Nutation  and  the  Secular 
Retardation  of  Earth’s  Rotation,  Phys.  Earth  Planet.  Interiors  1,  267-287, 
North-Holland  Publishing  Company,  Amsterdam,  1968. 

15.  P.  Musen  and  T.  Felsentreger,  On  the  Determination  of  the  Long  Period  Tidal 
Perturbations  in  the  Elements  of  Artificial  Earth  Satellites,  Celestial  Mechanics, 
Vol.  7,  No.  2,  pp.  256-279,  February  1973. 


32 


APPENDIX  A 

CONTRIBUTIONS  BY  TERMS  CONTAINING  THE  LOVE  NUMBER  hj 


CONTRIBUTIONS  BY  TERMS  CONTAINING  THE  LOVE  NUMBER  h. 


To  be  able  to  decide  which  terms  to  retain  and  which  ones  to  delete  from  the 
formula  for  the  tidal  surface  elevation,  the  following  (rather  simple  yet  sufficient) 
trial  case  was  done  by  hand  computation.  The  results  suggest  that  it  is  safe  at 
present  to  feature  in  the  computer  algorithm  only  the  two  terms  containing  the 
Love  number  hj  and  to  omit  all  the  others. 


g 

9.81 

m/sec^ 

Approximate  surface 
gravity 

G 

6.673  E - 08 

CTC?  jgx  sec^ 

Current  value 

R 

6378150 

m 

Typical  value 

4.9177  E+  18 

cm^  /sec^ 

Current  value 

Ms 

1.3291  E + 26 

cm^  /sec^ 

Current  value 

hj 

0.6 

Typical  values,  sources 

h3 

0.3 

. indicated  in  Appendix  C 

••l 

400  000 

km 

Typical  value 

■■s 

150  000  000 

km 

Typical  value 

•Hi.  >2  = Oi  f 

rl 

1 0"  ’ * P2  (cos  ) 

(201) 

= h,  7 

R^ 

rt 

1 0"  ' * Pj  (cos  y^^ ) 

(202) 

<Hs)2  - ^ 7 

R^ 

10"'*P2(cos  7j) 

(203) 

A-1 


(204) 


(Hs's  “ "J  f f lO-'^PjCcosTs) 


(Hl)2  = 0.191  m 

P2(cos7l) 

(205) 

(Hl)3  =0.001  m 

P3  (cos  7^ ) 

(206) 

(Hs)^  = 0.098  m 

P2(cos  75) 

(207) 

(Hs)3  = 2 . 10-<^ 

m P3  (cos  75 ) 

(208) 

First  Case 

Station  at  and  sun  and  moon  above  intersection  of  equator  and  Greenwich 
meridian. 

X = = 0 

= 7s  “ ® cos7j^  = cos7g  = 1 
-*•  Pj(cos7|^)  = Pj(cos75)  = 1,  i = 1,2 
= 0.191  m 
(Hl)3  =0.001  m 
(Hs)^  = 0.098  m 


I 


Second  Case 


Station  at  intersection  of  equator  and  Greenwich  meridian.  Moon  and  sun  at 
intersection  of  equator  and  180°  meridian. 


X = = 0 

= 7^  = 180°  -*■  cos7j^  = cos7g 


Pjf-D^+l,  P3(-1)--4 


= -1 


(Hl)2  = 0.191  m 

(Hl)3  = -0.006  m 


(H5)2  = 0.098  m 


(Hs)j  = 0.000  m 


H = 0.283  m 


I 


Third  Case 

Station  at  intersection  of  equator  and  Greenwich  meridian.  Moon  and  sun  at 
intersection  of  equator  and  +90°  meridian  or  -90°  meridian. 

7 = 1^  = 0 

7j^  = 7^  = ±90°  -*•  cos  7j^  = cos  7s  “ 0 

P2(0)  = -1/2,  PjfO)  = 0 

(Hl)2  = -0.096  m 

(H,)3  = 0 

(Hs)2  = -0.049  m 

(Hs)3  =_Q 

H » -0.145  111 
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APPENDIX  B 

A NOTE  CONCERNING  THE  APPLICATION  TO  CELEST 


A NOTE  CONCERNING  THE  APPLICATION  TO  CELEST 


During  the  early  computer  runs  made  while  developing  the  present  algorithm, 
the  tidal  surface  elevation  was  recomputed  for  each  data  line  on  the  CELEST 
observation  tape.  The  fact,  that  these  data  lines  follow  each  other  at  rather  short 
intervals,  makes  the  algorithm  a costly  augmentation  of  the  observation  tape.  On  the 
other  hand,  it  may  be  concluded  from  the  data  in  the  section  on  “QUANTITATIVE 
RESULTS”  that  it  is  unnecessary  to  update  H that  frequently  as  it  changes  quite 
slowly. 

The  CELEST  observation  tape  is  divided  into  data  blocks,  each  data  block 

representing  a satellite  pass.  In  actual  CELEST  applications,  two  types  of  data 
blocks  occur,  namely  those  that  are  about  20  min  long  and  others  which  are  several 
hours  in  length  (durations  in  excess  of  12  hr  being  possible).  The  former  are 

associated  with  passes  over  the  observation  stations  of  typical  geodetic  satellites.  The 
latter  arise  from  high  altitude  satellite  orbits. 

A simple  analysis  suggested  that  for  both  types  of  station  passes  it  is  equally 
suitable  to  calculate  H once  every  10  min.  For  the  short  passes  this  was  shown  by 
entering  the  trial  data  listed  in  Appendix  A into  the  equation  for  H 

(Equation  103).  The  observing  station  was  assumed  to  be  located  on  the  equator 

with  sun  and  moon  in  line  with  each  other  and  also  positioned  in  the  equator.  The 
tidal  elevation  was  then 

H(t)  = a(3  cos^  7“  1)  (301) 

a ^ 1 5 cm  . 


Assuming 


T « wt , (302) 


the  time  rate  of  H was 


dH  , ~ 

— « - 3acj  sin  2cjt  , 
dt 


(304) 


|H|  3a2i  2 cm/10  min 

' 'm  ax 


(305) 


For  a typical  long  duration  pass,  the  maximum  total  variation  of  H was 

3o(  45  cm  , - lift  • ,,ft^v 

« as  1.5  cm/ 10  min  , (306) 

6 hrs  6 hrs 

which  is  roughly  the  amount  estimated  for  the  20-min  pass.  Assuming  that  a 
discrepancy  of  2 cm  in  observing  station  elevation  is  yet  undetectable  by  the  means 
of  satellite  geodesy,  it  was  resolved  to  settle  on  the  10-min  interval  length  for  both 
types  of  satellite  passes. 
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NOTES  ON  THE  EFFECT  OF  THE  TIDAL  MASS  RELOCATION  ON 
TERRESTRIAL  GRAVITY 


The  present  analysis  of  the  tidal  surface  displacement  requires  that  the 
geometrical  Love  numbers  are  available.  In  particular,  h^  is  assumed  to  be  known  to 
about  ten  percent.  The  value  h2  = 0.6  was  specified  for  use  with  the  computer 
routine  for  the  surface  displacement.  It  was  adopted  from  Reference  11, 
Pages  27-29;  Reference  12,  Page  XVll;  Reference  13,  Page  368;  and  Reference  14, 
Page  285. 

For  the  purposes  of  the  computer  routine,  hj  is  irrelevant.  However,  to 
demonstrate  that  (which  is  done  in  Appendix  A),  a value  for  this  particular  Love 
number  had  to  be  known,  hj  = 0.3  was  finally  adopted  from  data  tabulated  in 
Reference  13,  Page  370. 

By  their  definition,^’*’''  the  geometrical  Love  numbers  and  the  dynamical  ones 
are  related  to  each  other.  For  this  reason,  it  was  desired  to  have  an  algorithm 
capable  of  deriving  a value  for  hj  from  the  Love  numbers  kj.  The  latter  are  already 
being  used  in  various  computer  programs*’'*  to  account  for  the  gravitational 
perturbation  which  the  tidal  bulge  of  the  solid  earth  exerts  on  satellites.  A suitable 
algorithm  was  obtained  as  follows. 

For  example,  as  mentioned  in  Reference  9,  Page  3,  the  tidal  deformation  of 
the  earth’s  surface  distorts  the  gravity  field.  Specifically,  the  presence  of  a particular 
solid  harmonic  in  the  tide  generating  potential  field  causes  a surface 

displacement 

H„  = h„  y (20) 

which  in  turn  produces  a distortion 


more  that  the  earth  is  spherical  (surface  radius  R)  and  that  its  density,  p,  is 
constant.  The  surface  area  element  is  dS.  The  surface  density  corresponding  to  the 
tidal  mass  redistribution  is  a. 


p HdS  = adS  , 


a = pH  , 


Considering  certain  relationships  from  the  main  body  of  this  report, 


f a^P^fsin  ^)  + 2 a p;;-  (sin  >fi)  cos  m X 

S K m - I 


^ m PIT  (si"  si" 

m=l  nm  n ' 


and  the  perturbing  potential  (acting  on  the  satellite),  T^^,  can  now  be  expressed  as  a 
Poisson  integral. 


T„  = G JJ 


r !>  a„(/,X')dS' 
(S)  IP-  r I 


(406) 


(T„)g  = (T„),=  r 


Once  more  using  Equations  (3b)  and  (6)  and  restricting  the  integration  to  computing 
surface  values  of  the  potential  only. 


(T„)s  = (Tn),  +(T„), 


~ J Jdi/j'dX'Yj,(/,X')cos^' ♦ 2 P^(sini^)P^(sini^') 
(S) 

26 

(T„)2  = J{dv»'dX'Y„(^',X')cos/ 

^ / C \ 


•»  V (v  “ u)l 

* -^0  J.  ^;;7^P:(sin^)P:;(sin/)cosMX-X') 


n - , P G 
" "i  R"-' 


(T„),  = 

4n 

2n  + 1 

2n 

n 

— B 

S a 

R " 

1 n)i 

PH  (sin  ifi)  cos  /iX 


■*■  J,  ’’nM  ^;;^PH(sinvJ)sinAiX 


and,  consequently, 


47r  p Y„(p,X)  . 


This  means  that  (Ti,)^  is  proportional  to  the  proportionality  factor  being  the 

Love  number  k , 

n 


<T»ls  = k„(U.),.,  . 


(415) 


F 


47r  p 

" 2n  + 1 g " 


!A  familiar  theorem  of  potential  theory  states  that  a harmonic  function  (potential)  is 

uniquely  determined  by  its  boundary  values  on  a specified  surface.  If  applied  to 

I Equation  (401),  the  relationship 

! 

I 

I'  T„(r,^,X)  = k„  U„(r,^,X)  (416) 


J 


I 


will  appear  justified.  This  latter  equation  is  the  one  which  is  frequently  invoked  to 
define  k^. 

If  one  enters  the  values  for  k^  and  kj  (actually  kj^  and  k^^,  or  averages. of 
the  variable  kj  and  k, ) listed  in  Reference  15,  Page  460  into  Equation  (415),  the 
following  values  result  for  hj  and  hj : 


kj  = 0.3  hj  = 0.9 

kj  = 0.1  hj  = 0.4 


(417) 
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While  the  two  values  for  the  dynamical  Love  numbers  are  generally  accepted  in 
magnitude,  the  two  h values  are  too  large.  This  failure  of  the  just  outlined  theory 
is  believed  to  be  caused  by  the  assumption  of  constant  density.  The  author  intends 
to  produce  a multilayer  model  which  he  hopes  will  yield  more  realistic  hj  values. 
Although  there  appears  at  present  no  need  to  refine  the  model  by  using  surface 
harmonic  expansions  for  the  hj  similar  to  those  used  in  Reference  11  for  the  kj,  it 
is  intended  to  try  that  while  deriving  the  improved  hj  = hj(kj)  algorithm. 
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